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ABSTRACT 

We report on Chandra and XMM-Newton observations of a sample of 13 quasars at z « 4.7-5.4 mostly taken 
from the Sloan Digital Sky Survey (SDSS). The present sample complements previous X-ray studies of z > 4 
quasars, in which the majority of the objects are optically more luminous and at lower redshifts. All but two of 
our quasars have been detected in the X-ray band, thus doubling the number of z > 4.8 X-ray detected quasars. 
The two non-detections are likely to be due to a short exposure time (SDSSp J033829. 31+002156. 3) and to the 
presence of intrinsic absorption (SDSSp J173744. 87+582829. 5). We confirm and extend to the highest redshifts 
the presence of a correlation between AB 1450(1 +z ) magnitude and soft X-ray flux for z > 4 quasars, and the presence 
of a steeper optical-to-X-ray spectral energy distribution (parameterized by a ox ) for high-luminosity, high-redshift 
quasars than for lower-luminosity, lower-redshift quasars. The second effect is likely due to the known anti- 
correlation between oc ox and rest-frame 2500 A luminosity, whose significance is confirmed via partial correlation 
analysis. The joint w 2.5-36 keV rest-frame spectrum of the z > 4.8 SDSS quasars observed thus far by Chandra 
is well parameterized by a power-law with photon index T= 1 .84^q lh this photon index is consistent with those of 
z ~ 0-3 quasars and that obtained from joint spectral fitting of z ~ 4.1-4.5 optically luminous Palomar Digital Sky 
Survey quasars. No evidence for widespread intrinsic X-ray absorption has been found (Na ~ 4.0 x 10 22 cm~ 2 
on average at 90% confidence). We also obtained Hobby-Eberly Telescope (HET) photometric observations for 
eight of our target quasars. None of these shows significant (> 30%) optical variability over the time interval of a 
few years (in the observed frame) between the SDSS and HET observations. 



1. INTRODUCTION 

One of the main themes in modern astronomy is the study 
of the earliest massive objects to form in the Universe, now 
known as far back as the end of the reionization epoch (e.g., 
Rees 1999; Loeb & Barkana 2001; Fan et al. 2002; Hu et al. 
2002). Quasars at z 4 can be studied in detail; because of their 
large broad-band luminosities, they are typically detectable at 
most wavelengths [e.g., in the radio, Schmidt et al. 1995; Stern 
et al. 2000; in the millimeter, Omont et al. 2001; Carilli et al. 
2001; in the sub-millimeter, McMahon et al. 1999; Isaak et al. 
2002; in the near-infrared, Barkhouse & Hall 2001; Pentericci 
et al. 2003; in the optical, Schneider, Schmidt, & Gunn 1989; 
Fan et al. 2001; in the X-ray band, Kaspi, Brandt, & Schneider 
2000 (hereafter KBS00); Vignali et al. 2001, 2003 (hereafter 
VOla, V03); Bechtold et al. 2003]. Ground-based, large-area 
optical surveys such as the Sloan Digital Sky Survey (SDSS; 
York et al. 2000), the Automatic Plate Measuring facility sur- 
vey (BRI; Irwin, McMahon, & Hazard 1991), and the Palomar 
Digital Sky Survey (PSS; Djorgovski et al. 1998) have discov- 
ered more than 300 quasars at z > 4 thus far up to z = 6.43 (Fan 
et al. 2003). 6 

To extend our knowledge of the X-ray properties of z > 4 
quasars, we have started a program to observe, with Chan- 
dra and XMM-Newton, both the optically most luminous 
(M B « -28.4 to -30.2) z~ 4.1-4.5 PSS/BRI quasars (e.g., 
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V03) and the higher-redshift, optically fainter SDSS quasars 
(VOla; Brandt et al. 2002; see Brandt et al. 2003 for a recent 
review). Selecting z > 4 quasar samples for X-ray observations 
in two different regions of the luminosity-redshift parameter 
space allows study of any dependence of the X-ray properties 
upon the optical selection criteria. In particular, it is possible to 
test for dependencies of z > 4 quasar X-ray continua on optical 
luminosity or redshift. 

Here we report the results obtained from 1 1 Chandra obser- 
vations and one XMM-Newton observation of a sample of 13 
high-redshift (z « 4.7-5.4), moderately luminous (Mb ~ —26.1 
to —28.5) quasars, mostly discovered by the SDSS and pub- 
lished in Anderson et al. (2001). One of these quasars (RD 657; 
Djorgovski et al. 2003) was serendipitously observed and de- 
tected in the XMM-Newton observation (i.e., there were two 
objects in one field). The only target quasar not coming from 
the SDSS data is BR B0305-4957 (Storrie-Lombardi et al. 
2001; Peroux et al. 2001). Given that the accretion rate 
(relative to the Eddington rate) of quasars may increase with 
redshift (e.g., Kauffmann & Haehnelt 2000), the first quasars 
could show unusual phenomena such as "trapping radius" ef- 
fects (e.g., Begelman 1979) or accretion-disk instabilities (e.g., 
Lightman & Eardley 1974). Furthermore, our targets are bright 
enough to assure in most cases X-ray detections with snap- 
shot (e.g., « 4-10 ks) observations, given the correlation found 
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' See http://www.astro.caltech.edu/~george/z4.qsos for a listing of known high-redshift quasars. 
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between AB 1450 ( 1+?; ) magnitude and soft X-ray flux (see, e.g., 
VOla; V03). The number of z > 4.8 X-ray detected quasars is 
approximately doubled by the objects presented here. As shown 
in Figure 1, where the AB 1450 ( 1+Z \ magnitude is plotted as a 
function of redshift, the objects presented in this paper suitably 
complement previous X-ray studies of z > 4 quasars. All but 
one of these quasars are likely radio-quiet (see §3 for details). 




4 4.5 5 5.5 6 
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FIG. 1. — AB 1450 (i +z ) versus redshift for the known z > 4 AGNs [small 
dots; taken from Djorgovski as of 2002 July (see Footnote 6)]. The X-ray de- 
tected (undetected) quasars presented in this paper are shown as filled (open) 
triangles (Chandra observations) and a filled (open) square (XMM-Newton 
observation). Filled and open circles indicate z > 4 X-ray detected AGNs and 
tight upper limits, respectively, from previous observations [Schneider et al. 
1998; Waddington et al. 1999 (note that this identification is still debated; 
L. L. Cowie 2002, private communication); KBS00; VOla; Brandt et al. 2001, 
2002; Silverman et al. 2002; Vignali et al. 2002 (V02); Barger et al. 2002; 
Bechtold et al. 2003; V03; Castander et al. 2003], most of which were per- 
formed by Chandra and ROSAT. 

Hereafter we adopt Hq=70 km s -1 Mpc -1 in a A-cosmology 
with i2 M =0.3 and il A =0.7 (e.g., Turner 2002). 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Chandra Observations 

Eleven of the 13 quasars were observed by Chandra dur- 
ing Cycle 3 with « 3.7-9.8 ks observations (see Table 1; here- 
after all of the quasars will be referred to via their abbreviated 
names). All of the sources were observed with the Advanced 
CCD Imaging Spectrometer (ACIS; Garmire et al. 2003) with 
the S3 CCD at the aimpoint. Faint mode was used for the 
event telemetry format, and ASCA grade 0, 2, 3, 4 and 6 events 
were used in the analysis. This particular grade selection is 
standard and appears to optimize the signal-to-background ra- 
tio (see §6.3 of the Chandra Proposers' Observatory Guide). 
Source detection was carried out with WAVDETECT (Freeman 
et al. 2002). We used wavelet transforms (using a Mexican hat 
kernel) with wavelet scale sizes of 1, 1.4, 2, 2.8, and 4 pixels 
and a false-positive probability threshold of 10~ 4 . Given the 
small number of pixels being searched due to the known source 

7 To account for the quantum efficiency decay of ACIS at low energies, possibly 
dependent correction to the ACIS quantum efficiency. We generated new ancilk 
throughout this paper. 



positions and the excellent angular resolution of Chandra, the 
probability of spurious detections is extremely low; most of 
the sources were in fact detected at a false-positive probabil- 
ity threshold of 10~ 6 . Typically, detections were achieved with 
a wavelet scale of 1.4 pixels or less. 

Source searching was performed in four energy ranges for 
consistency with our previous work (VOla; V03): the ultra- 
soft band (0.3-0.5 keV), the soft band (0.5-2 keV), the hard 
band (2-8 keV), and the full band (0.5-8 keV); in the redshift 
range of z ~ 4.7-5.4 covered by our sample, the full band cor- 
responds to the w 2.8-51 keV rest-frame band. All but one of 
the sources observed by Chandra were detected. The Chan- 
dra positions of the detected quasars were found to lie within 
0.2"-0.8" of their optical positions (see Table 1); this is consis- 
tent with the expected positional error. Given the X-ray weak- 
ness of SDSS 0756+4502 (two full-band counts), we will dis- 
cuss the reliability of its detection in more detail in §2.1.1. The 
Chandra non-detection, SDSS 1737+5828, will be discussed 
further in §4. 1 . 

The WAVDETECT photometric measurements are shown in 
Table 2; they are consistent with independent measurements 
obtained using manual aperture photometry and a 2" radius cir- 
cular cell. Background is typically negligible in the source ex- 
traction region, except for SDSS 0941+5947, where the obser- 
vation is partially affected by a background flare; this is taken 
into account in Table 2. The only object with enough counts 
to derive a useful power-law photon index via direct X-ray 
spectral analysis is SDSS 0231—0728. Assuming only Galac- 
tic absorption, it has Y = 2.4+Qg using the unbinned data and 
the C-statistic (Cash 1979) with xspec Version 1 1 .2.0 (Arnaud 
1996); errors are at the 90% confidence level for one interest- 
ing parameter (AC = 2.71; Avni 1976; Cash 1979). 7 This pho- 
ton index is consistent, within the errors, with those of z ~ 0-3 
quasars in the rest-frame 2-10 keV band (r « 1.7-2.3; e.g., 
George et al. 2000; Mineo et al. 2000; Reeves & Turner 2000; 
Page et al. 2003). 

The soft-band images of the detected quasars (full-band im- 
age for BR 0305-4957), adaptively smoothed using the algo- 
rithm of Ebeling, White, & Rangarajan (2003), are shown in 
Figure 2. 

2.1.1. Chandra Detection of SDSS 0756+4502 

Similarly to the z = 5.27 quasar SDSS 1208+0010 presented 
in §2.2 of VOla, SDSS 0756+4502 (z = 4.80) shows only two 
X-ray photons (both in the observed 0.5-2 keV band) within 
0.7" of the optical position of the source. The source is de- 
tected by WAVDETECT using a false-positive threshold of 10~ 4 . 
To investigate further whether the SDSS 0756+4502 X-ray de- 
tection can be considered real, we extracted a 400x400 pixel 2 
(« 200x200 arcsec 2 ) region centered on the optical position 
of the quasar, excluding the immediate vicinity of the quasar 
itself. We covered this region with 30,000 circles of 1" ra- 
dius whose centers were randomly chosen. The choice of the 
1" radius aperture appears reasonable, since the full-band 80% 
encircled-energy radius for on-axis sources is w 0.69" (see 
Fig. 6.3 of the Chandra Proposers' Observatory Guide) and 
the average positional error for the other quasars detected by 
Chandra is w 0.5". The counts obtained in each circle were 
averaged, giving 0.060 counts circle -1 in the full band and 

caused by molecular contamination of the ACIS filters, we have applied a time- 
ry response files which have been used in the X-ray spectral analyses presented 
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TABLE 1 

Chandra AND XMM-Newton OBSERVATION LOG 



Object 




Optical 


Optical 


Aopt-x d 


X-ray 


Exp. Time " 


HET 




Name 


z 


«2000 


§2000 


(arcsec) 


Obs. Date 


(ks) 


Obs. Date 


Ref. 


SDSSp J020651.37+121624.4 


4.81 


02 06 51.4 


+ 12 16 24.4 


0.8 


2001 Dec 29 


5.88 


2001 Dec 10 


(1) 


SDSSp J023137.65-072854.5 


5.41 


02 31 37.6 


-07 28 54.5 


0.2 


2002 Sep 27 


4.20 




(1) 


BRB0305-4957 


4.73 


03 07 22.9 


-49 45 47.8 


0.4 


2002 Sep 22 


3.71 




(2,3) 


SDSSp J033829. 31+002156. 3 c 


5.00 


03 38 29.3 


+00 21 56.3 




2002 Feb 22 


5.49 




(4) 


RD 657 c 


4.96 


03 38 30.0 


+00 18 39.9 


1.7 


2002 Feb 22 


4.85 




(5) 


SDSSp J075618.14+410408.6 


5.09 


07 56 18.1 


+41 04 08.6 


0.2 


2002 Feb 08 


7.33 


2001 Dec 09 


(1) 


SDSSp J075652.07+450258.9 


4.80 


07 56 52.1 


+45 02 58.9 


0.5 


2002 Oct 04 


6.97 


2001 Dec 10 


(1) 


SDSSp J091316.56+591921.5 


5.11 


09 13 16.6 


+59 19 21.5 


0.7 


2002 Mar 07 


9.84 


2002 Feb 07 


(1) 


SDSSp J094108.36+594725.8 


4.82 


09 41 08.4 


+59 47 25.8 


0.2 


2001 Oct 25 


4.16 


2001 Nov 16-18 


(1) 


SDSSp J095 15 1.17+594556.2 


4.86 


09 51 51.2 


+59 45 56.2 


0.3 


2002 Apr 14 


5.13 


2002 Feb 07 


(1) 


SDSSp J102332.08+633508.1 


4.88 


10 23 32.1 


+63 35 08.1 


0.1 


2002 Oct 09 


4.70 




(1) 


SDSSp J173744.87+582829.5 


4.94 


17 37 44.9 


+58 28 29.5 




2002 Aug 05 


4.62 


2002 Jun 08 


(1) 


SDSSp J221644.02+001348.3 


4.99 


22 16 44.0 


+00 13 48.3 


0.4 


2002 Nov 01 


7.41 


2002 Jun 14 


(1) 



NOTE. — The optical information for the quasars presented here can be found in the papers cited in the reference column. For BR B0305— 4957 we 
use the redshift z = 4.73 from Peroux et al. (2001) derived from two emission lines (Si IV-O IV] and C IV) in their high-quality optical spectrum. 
Storrie-Lombardi et al. (2001) reported z = 4.78 (see their Table 1) and z = 4.82 (see their Fig. 3) derived from a lower-quality optical spectrum 
with less wavelength coverage. 

a Distance between the optical and X-ray positions; a blank entry indicates no X-ray detection. b The Chandra exposure time has been corrected for 
detector dead time. The XMM-Newton exposure time, computed from the exposure maps (i.e., taking into account the decreasing effective exposure 
at increasing off-axis angles from the aimpoint), is for the observation after the background flares have been removed (see §2.2). c XMM-Newton 
observation. 

REFERENCES. — Anderson et al. 2001; (2) Peroux et al. 2001; (3) Storrie-Lombardi et al. 2001; (4) Fan et al. 1999; (5) Djorgovski et al. 2003. 



0.022 counts circle -1 in the soft band. The Poisson probabil- 
ities of obtaining two counts or more when 0.060 and 0.022 
counts are expected are w 1.7 x 10~ 3 and w 2.4 x 10~ 4 , cor- 
responding to one-tailed Gaussian probabilities of w 3.3c and 
w 3.9a, respectively. 

We have directly verified the validity of the Poisson ap- 
proximation for the SDSS 0756+4502 field using the count 
distribution from the 30,000-circle analysis described above. 
The results of our simulations are in reasonable agreement 
with the Poisson probabilities. Therefore, in the following, 
SDSS 0756+4502 will be considered detected by Chandra. 



is found in each quasar field. A comparison of their 0.5-2 keV 
cumulative number counts (excluding, of course, the quasars 
themselves) with those obtained from the ROSAT deep sur- 
vey of the Lockman Hole (Hasinger et al. 1998) and the 
XMM-Newton moderately deep survey HELLAS2XMM (Baldi 
et al. 2002) indicates that in the proximity of the SDSS 
quasars there is no evidence for an X-ray source overdensity. 
We found a surface density of N(> S)=751+3q[ deg~ 2 for the 
X-ray sources surrounding our SDSS targets (where S is the 
soft X-ray flux limit we reached with the present observations: 
«2 x 10~ 15 ergcirr 2 s" 1 ). 



2.1.2. X-ray Sources Close to the Quasars 

Following the method used in §3.4 of V03, we have 
searched for possible companions for all of the objects ob- 
served by Chandra over a region of « 100"xl00" centered 
on the quasars (see Fig. 2), corresponding to a physical scale 
of w 630 x 630 kpc 2 at the average redshift of our sample 
(z « 5.0). Approximately one faint serendipitous X-ray source 



2.1.3. X-ray Extension of the Quasars 

To constrain the presence of gravitational lensing (e.g., 
Barkana & Loeb 2000; Comerford, Haiman, & Schaye 2002; 
Wyithe & Loeb 2002), X-ray jets close to the quasars (e.g., 
Schwartz 2002), and X-ray scattering by dust in the intergalac- 
tic medium (e.g., Telis et al. 2003), we performed an analy- 
sis of X-ray extension for all the quasars observed by Chandra 



Table 1 

Chandra SOURCES: X-RAY COUNTS 







X-ray Counts" 




Object 


[0.3-0.5 keV] 


[0.5-2 keV] 


[2-8 keV] 


[0.5-8 keV] 


SDSS 0206+1216 


<3.0 




<4.8 


6 9+ 18 


SDSS 0231-0728 


2.0 + _\l 


22.0+58 


2Xt\l 


23.9+|° 


BR 0305-4957 


<3.0 


2.0+fJ 


2.0t\l 


4.0+t 2 . 


SDSS 0756+4104 


2 0+ 2 ' 7 

Z.U_j 3 


10.0+ 4 / 3 


3.9I 3 ; 2 


14.8+ 4 ^ 


SDSS 0756+4502 


<4.8 


2.0+!? 


< 3.0 


2.0+fl 


SDSS 0913+5919 


<4.8 




<4.8 


6 0+ 36 

O.U_ 24 


SDSS 0941+5947 


<3.0 


4.8+J;? 


<3.0 


A o+ 3 ' 2 
^■ z -2.0 


SDSS 0951+5945 


<3.0 




<4.8 


-t q+3.2 


SDSS 1023+6335 


<3.0 


3 0+ 2 ' 9 

■'•"-1.6 


<4.8 


4 0+ 3 ' 2 


SDSS 1737+5828 


<3.0 


<3.0 


<3.0 


<3.0 


SDSS 2216+0013 


<4.8 


6 0+ 3 6 

O.U_ 24 


<4.8 


6 n+3.8 



" Errors on the X-ray counts were computed according to Tables 1 and 2 of Gehrels 
(1986) and correspond to the la level; these were calculated using Poisson statistics. 
The upper limits are at the 95% confidence level and were computed according to 
Kraft, Burrows, & Nousek (1991). 
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FIG. 2.— Chandra 0.5-2 keV (0.5-8 keV for BR 0305-4957) and XMM-Newton 0.3-12 keV (for RD 657) images of the X-ray detected quasars presented in 
this paper. In each panel, the horizontal axis shows the Right Ascension, and the vertical axis shows the Declination (both in J2000 coordinates). North is up, and 
East to the left. The images are 97" x 97" and have been adaptively smoothed at the 2a level [except for BR 0305-4957, RD 657, SDSS 0951+5945 (1.5o level), 
and SDSS 0756+4502 (1.3a level)] using the algorithm of Ebeling et al. (2003). Crosses mark the optical positions of the quasars. The image of RD 657 does not 
cover the position of SDSS 0338+0021, since the two quasars are separated by 196". 



following the procedure described in §3.5 of V03. All of the 
quasars' X-ray images are consistent with the sources being 



point-like. The separation between multiply imaged sources 
at z « 4-5 is expected to be w 1-2" (Barkana & Loeb 2000). 
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Typical flux ratios of multiply imaged quasars at high redshift 
are « 2-8 (G. Chartas 2003, private communication), being 
strongly dependent on the shape, mass, redshift, and location of 
the lensing system. Given the limited number of counts for our 
detected sources, we expect that only in w 30% of the present 
sample would we be able to detect such multiple images. 

2.1.4. X-ray Variability of the Quasars 

Given tentative reports of increasing quasar X-ray variabil- 
ity with redshift (in the sense that quasars of the same X-ray 
luminosity are more variable at z > 2; Manners, Almaini, & 
Lawrence 2002), we have searched for X-ray variability by ana- 
lyzing the photon arrival times in the 0.5-8 keV band using the 
Kolmogorov-Smirnov test. No significant variability was de- 
tected from any of our quasars. We note, however, that the anal- 
ysis of X-ray variability for the quasars studied here is limited 
by both the short X-ray exposures in the rest frame (« 15 min; 
see Table 1) and the low X-ray fluxes of most of our quasars. 

2.2. XMM-Newton Observation 

The SDSS 0338+0021 field was observed by XMM-Newton 
for a nominal exposure time of 20.0 ks with the European Pho- 
ton Imaging Camera (EPIC) pn (Striider et al. 2001) and for 
22.3 ks with the EPIC MOS (Turner et al. 2001) instruments. 
The data obtained with the Reflection Grating Spectrometer 
(RGS; den Herder et al. 2001), with a nominal exposure time of 
22.8 ks, will not be considered in the following given the weak- 
ness of the target. Both the pn and MOS data were acquired 
in the "full-frame mode" using the thin filter (due to the likely 
weakness of the target in the X-ray band). The data were re- 
duced with the XMM-Newton Science Analysis Software (SAS, 
Version 5.3.3) with the latest calibration products; we used the 
EPCHAIN and EMPROC tasks to generate valid photon lists. We 
filtered the data to avoid background flares (due tow 100 keV 
protons; e.g., De Luca & Molendi 2002) which affect most of 
this XMM-Newton observation. To select good time intervals, 
we apply thresholds of 0.35 counts s _1 and 0.15 counts s _1 to 
pn and MOS data in the 10-12.4 keV and 10-13 keV energy 
bands, respectively (A. Baldi 2002, private communication). 
In these energy bands, the counts are mostly from the back- 
ground. The final, on-axis useful exposure time for this obser- 
vation is 5.5 ks forpn and MOS. Both single-pixel and double- 
pixel events (patterns 0-4) were used when extracting the pn 
counts, while patterns 0-12 were used for the MOS. 

We performed source detection in the 0.3-2 keV, 0.5-2 keV, 
0.5-10 keV, 2-10 keV, and 0.3-12 keV bands using the tasks 
EBOXDETECT and EMLDETECT and adopting a detection like- 
lihood threshold of 10, corresponding to the « 4a detection 
level. We did not consider the energy range below 0.3 keV be- 
cause of the enhanced background at low photon energies. The 
target, SDSS 0338+0021, was detected in none of the bands; 
its X-ray non-detection will be discussed further in §4. 1 . 

After our XMM-Newton observation was assigned, Djorgov- 
ski et al. (2003) found a z = 4.96 quasar (RD 657) located 
196" from the optical position of SDSS 0338+0021. In the 
XMM-Newton pn field-of-view we found a 4.3a X-ray source 
detected 1.7" from the optical position of RD 657 (see Table 1) 



in the 0.3-12 keV band (15.3 counts). We note the XMM- 
Newton positional error is larger (« 4-5"; Metcalfe 2002) than 
that of Chandra; therefore the X-ray and optical positions are 
consistent. Assuming the integral 0.5-2 keV number counts 
of Hasinger et al. (1998) and an X-ray error circle of 1.7", 
we expect « 2.5 x 10~ 4 spurious associations given the soft 
X-ray flux of this source (3.6xl0~ 15 erg cm~ 2 s _1 ; see Ta- 
ble 3). Therefore, we are confident that the serendipitously 
detected X-ray emission is associated with RD 657. For this 
source we compared the 4.3a detection obtained by the SAS 
detection tasks with that obtained in a different way with S EX- 
TRACTOR (Bertin & Arnouts 1996), finding similar results in 
both the number of counts and the significance of the source 
(3.6a). 

Inspection of the X-ray image at the source position con- 
firms the X-ray detection (see Fig. 2 for the 1.5a level adap- 
tively smoothed image in the 0.3-12 keV band). RD 657 also 
seems to be present in the 0.3-2 keV adaptively smoothed im- 
age, although the quasar is not detected in this band by the 
SAS procedure described above. Even though there has been 
speculation that SDSS 0338+0021 and RD 657 mark a large- 
scale structure, perhaps a proto-cluster at z ~ 5 (Djorgovski 
et al. 2003), we found no clear evidence for diffuse X-ray 
emission or an over-density of X-ray sources around these two 
quasars. Our constraints are not tight, however, due to the 
short exposure time and relatively high background of XMM- 
Newton. Since the background in the XMM-Newton obser- 
vation is higher and more complex than in the Chandra ob- 
servations, for SDSS 0338+0021 and RD 657 we do not re- 
port the counts in individual bands but directly derived the 
main X-ray parameters in Table 3 (see §3). We have also 
analyzed the XMM-Newton data with less conservative back- 
ground screening choices and our main results are unchanged; 
SDSS 0338+0021 is not detected and RD 657 is detected. 

2.3. Hobby-Eberly Telescope Observations 

Given that quasars frequently vary in both the optical and 
X-ray bands, one must be concerned about the reliability 
of the frequently employed optical-to-X-ray power-law slope 
whenever observations in the two bands are not simultaneous. 
To search for optical variability, eight of the 13 quasars in 
our sample were observed with the 9-m Hobby-Eberly Tele- 
scope (Ramsey et al. 1998). 8 Short (2 min; 4 min for 
SDSS 0941+5947) /-band images of the quasar fields were 
taken with the HET's Marcario Low-Resolution Spectrograph 
(LRS; Hill et al. 1998a,b; Cobos Duenas et al. 1998) generally 
within 1-2 months of the X-ray observations (1-2 weeks in the 
rest frame; see Table l). 9 The HET images covered fields « 4' 
on a side, and using the published finding charts of the fields it 
was possible to obtain an approximate photometric calibration 
using stellar objects in the fields. Although it is difficult to give 
precise estimates given the difference between the LRS and 
SDSS bandpasses, it is clear that none of the quasars observed 
with the HET has varied significantly (by £ 30%) in brightness 
from its published value. These constraints on the optical flux 
variability allow for at most an uncertainty of « 0.05 in the a ox 
determination (see §3). The lack of large brightness variations 
in the sample is not surprising considering the relatively small 



8 HET observations were not performed for SDSS 0231-0728, BR 0305-4957, SDSS 0338+0021, RD 657, and SDSS 1023+6335. 

9 The only exception is SDSS 0756+4502, observed by the HET 10 months before the X-ray observation. 
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amount of time in the rest frame (several months) between the 
SDSS photometry and the X-ray/HET observations. 

To investigate the nature of SDSS 1737+5828, the only 
quasar that was not detected in the Chandra observations, we 
performed a 30-min LRS spectroscopic observation. This was 
taken with the 300 line mm -1 grating, an OG515 blocking fil- 
ter, and a slit width of 2" (with a seeing of w 1.6"), providing 
a spectrum in the wavelength range 5100-10200 A at a resolv- 
ing power of w 240-300. The spectroscopic data were reduced 
using standard IRAF routines. The wavelength calibration was 
performed using a Ne arc lamp, while the flux calibration was 
achieved using the standard star BD +17 4708, observed dur- 
ing the same night. The HET spectrum (black curve in Fig. 3) is 
discussed in §4.1 along with the SDSS spectrum (adapted from 
Anderson et al. 2001; grey curve in Fig. 3). 




7000 8000 9000 

Wavelength (A) 

FIG. 3.— Optical spectra of SDSS 1737+5828 obtained with the HET (30- 
min exposure; black curve) and by the SDSS (adapted from Anderson et al. 
2001; grey curve). No correction for the telluric lines has been applied to the 
HET spectrum (the affected regions of the HET spectrum are shaded). The 
likely absorption features from N V and Lya are shown (see §4.1 for details). 

A direct comparison of the continuum levels of the two spec- 
tra is difficult. While the uncertainty in the calibration of the 
SDSS spectrum is less than w 20%, it is hard to obtain an ac- 
curate HET flux calibration mostly because of the variations of 
effective aperture produced by its optical design (in addition to 
problems due to slit losses). 

3. MULTI-WAVELENGTH PROPERTIES OF Z > 4.7 QUASARS 

The principal optical, X-ray, and radio properties of the 
quasars in our sample are shown in Table 3. A description is 
as follows: 

Column (1 ). — The abbreviated name of the source. 

Column (2). — The Galactic column density (from Dickey & 

Lockman 1990) in units of 10 20 cm~ 2 . 

Column (3). — The monochromatic rest-frame AZ? 1450 ( 1+Z ) 
magnitude (defined in §3b of Schneider et al. 1989) with es- 
timated errors m 0.1 mag. The absolute photometry is correct 
to within w 0.05 mag (Anderson et al. 2001), while the compu- 
tation of the Afi 14 5o( 1+z ) magnitude using the photometric mea- 
surements and the composite SDSS quasar spectrum (Vanden 



Berk et al. 2001) is correct to within w 0.05-0.10 mag. 
Columns (4) and (5). — The 2500 A rest-frame flux density 
and luminosity. These were computed from the AB 1450 ( 1+Z ) 
magnitude assuming an optical power-law slope of a = —0.79 
(S v <x v a ; Fan et al. 2001) to allow direct comparison with the 
results presented in VOla and V03. The 2500 A rest-frame 
flux densities and luminosities are reduced by w 15% when 
the power-law slope of the optical continuum is changed to 
a = -0.5 (e.g., Vanden Berk etal. 2001; Schneider etal. 2001). 
Column (6). — The absolute B-band magnitude computed us- 
ing a = —0.79. If a = —0.5 is adopted for the extrapolation, 
the absolute B-band magnitudes are fainter by « 0.35 mag. 
Columns (7) and (8). — The observed count rate in the 
0.5-2 keV band and the corresponding flux (Fx), corrected 
for Galactic absorption. This flux has been calculated using 
PIMMS (Version 3. 2d; Mukai 2002) and a power-law model 
with r = 2.0, as derived for samples of z ~ 0-3 quasars (e.g., 
George et al. 2000; Mineo et al. 2000; Reeves & Turner 2000; 
Page et al. 2003) and z > 4 optically luminous, PSS quasars 
(V03); see also §4.3. 

The X-ray fluxes reported in this paper have been corrected 
for the ACIS quantum efficiency decay at low energy (see 
Footnote 7). For each quasar, we produced an ancillary re- 
sponse file (ARF) assuming an exponential decay for the ACIS 
quantum efficiency. From each ARF we derived an effec- 
tive area that was used to replace the original ACIS-S effec- 
tive area in PIMMS. The 0.5-2 keV flux correction is w 20% 
for the Chandra observations presented here. All of the pa- 
rameters derived from the measured soft X-ray flux have been 
corrected. Following this procedure, we also corrected all 
of the X-ray fluxes of the z > 4 quasars previously observed 
by Chandra (see Table Al in the Appendix for an updated 
version of the relevant numbers for z > 4 quasars; also see 
http://www.astro.psu.edu/users/niel/papers/). 

The soft X-ray flux changes by « 1-2% over a reasonable 
range in the adopted power-law photon index, T = 1 .7-2.3. The 
soft X-ray flux derived from the 0.5-2 keV counts is generally 
similar (to within w 5-20%) to that derived from the 0.5-8 keV 
counts. 

For RD 657, the 0.5-2 keV flux and subsequent X-ray pa- 
rameters reported below have been derived using the counts in 
the 0.3-12 keV band, where it was detected. 
Columns (9) and (10). — The rest-frame 2 keV flux density 
and luminosity, computed assuming F = 2.0. 
Column (11). — The 2-10 keV rest-frame luminosity, cor- 
rected for Galactic absorption. 

Column (12). — The optical-to-X-ray power-law slope, a ox , 
defined as 

log(/2kev// 2 500A) 



a 



log(V2kev/V 2 500A) 



(1) 



where fa k e v and / 2 500 A are tne rest-frame flux densities at 
2 keV and 2500 A, respectively. The « la errors on a ox have 
been computed following the "numerical method" described in 
§ 1.7.3 of Lyons (1991). Both the statistical uncertainties on 
the X-ray count rates and the effects of possible changes in the 
X-ray (r w 1.7-2.3) and optical (a w -0.5 to -0.9; Schneider 
et al. 2001) continuum shapes have been taken into account (see 
§3 of VOla for further details). Changing the power-law slope 
of the optical continuum from a = —0.79 to a = —0.5 induces 
a small increase in the a ox values (more positive by w 0.03). 
Column (13). — The radio-loudness parameter (e.g., Keller- 
mann et al. 1989), defined as R=fs GHz/ /4400 A ( rest frame). 



Table 2 

Optical, X-ray, and Radio Properties of z > 4.7 Quasars Observed by Chandra and XMM-Newton 



Object 
(1) 


(2) 


(3) 


^2500^ 

(4) 


log(vLv)2500 

(5) 


M B 
(6) 


Count rate 
(V) 


/x" 

(8) 


fl keV 

(9) 


log(vLv)2 keV 

(10) 


log(L2-10 kev) f 
(11) 


OV 8 

(12) 


(13) 


SDSS 0206+1216 


6.82 


19.7 


0.74 


46.6 


-27.5 


i s\ s\ -4-0 f, 1 

i.oo + g;j 


4A -fl 


3.83 


44.9 


45.1 


-1.64 + 0.10 


< 19.2' 


SDSS 0231-0728 


3.05 


19.3 


1.06 


46.8 


-28.1 


5.24+{;g 




20.36 


45.7 


45.9 




< 4.2' 


BR 0305—4957 


2.05 


18.6 


2.03 


47.0 


-28.5 


0.54+g™ 


2 1 +29 
z ' 1 -1.4 


1.82 


44.5 


44.8 


— 1 Q4+ 016 


"' k . 


SDSS 0338+0021 


8.11 


20.0 


0.56 


46.5 


-27.2 


< 1.88 


< 3.2 


< 2.89 


< 44.8 


< 45.0 


< -1.65 


< 26.0' 


RD 657 


8.05 


21.1 


0.20 


46.0 


-26.1 


< 2.80 


3 6+ 1 ' 21 


3.21 


44.8 


45.0 


-1.46 + 0.08 


<71.3 i 


SDSS 0756+4104 


4.76 


20.0 


0.56 


46.5 


-27.3 


1 36+ - 59 

i - 5D -0.42 


5 7+ 2 - 5 

■'■'-1.8 


5.18 


45.1 


45.3 


-1.55 + 0.09 


< 7.3' 


SDSS 0756+4502 


4.95 


20.0 


0.56 


46.4 


-27.2 


29+ - 35 


1 2 +L7 

'•^-0.8 


1.07 


44.3 


44.5 


i ki+ 016 

-1O1-0.18 


< 6.3' 


SDSS 0913+5919 


3.81 


20.3 


0.42 


46.4 


-27.0 


S1+ - 34 
UJ 1 -0.29 


2 1 +L4 


1.89 


44.6 


44.8 


-1.67 + 0.11 


410.8 s 


SDSS 0941+5947 


2.16 


19.2 


1.17 


46.8 


-28.0 


, ,,+0.82 

1 - 1J -o.50 


4 4+32 


3.84 


44.9 


45.1 


-1.72 + 0.11 


< 12.2 1 


SDSS 0951+5945 


1.44 


19.6 


0.81 


46.6 


-27.6 


S9+ - 57 


2.2± 2 ;| 


1.96 


44.6 


44.8 


i -77+0.13 
— 1-'/_0.14 


< 17. t 


SDSS 1023+6335 


1.01 


19.5 


0.88 


46.7 


-27.7 


64+ 62 


2.41- 


2.15 


44.6 


44.8 


i 77 +0.13 
'■"-0.14 


< 16.? 


SDSS 1737+5828 


3.58 


19.1 


1.28 


46.8 


-28.1 


<0.65 


<2.7 


<2.37 


<44.7 


<44.9 


< -1.82 


< 3.0' 


SDSS 2216+0013 


4.68 


20.1 


0.51 


46.4 


-27.1 


0.81jjg 


3 5+ 20 


3.09 


44.8 


45.0 


-1.62 + 0.10 


< 7.7' 





Note. — Luminosities are computed using Hq-70 km s _1 Mpc -1 , £1^=0.3, and £2^=0.7. 

a From Dickey & Lockman (1990) in units of 10 20 cm -2 . 

b Resi-frame 2500 A flux density in units of 10 -27 erg cm -2 s~ 1 Hz~ 1 . 

c Observed count rale computed in the 0.5-2 keV band in units of 10 -3 counts s _1 . 

^Galactic absorption- corrected flux in the observed 0.5-2 keV band in units of 10 — erg cm -2 s — ' . These fluxes and the following X-ray parameters have been corrected for the ACIS quantum efficiency decay at low energy. 
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1 Rest-frame 2-10 keV luminosity corrected for Galaclic ahsorplion in units of erg s 1 . 

§ Errors have been computed following the "numerical method" described in § 1.7.3 of Lyons (1991); both the statistical uncertainties on the X-ray count rales and the effects of the observed ranges of Ihe X-ray and optical continuum shapes have been taken into 
account (see the text for details; see also §3 in VOla). 

'''4400^ 

' 1 .4 GHz flux density from NVSS (Condon et al. 1998). 
J 1 .4 GHz flux density from FIRST (Becker el al. 1 995). 
k No radio data are available. 

'The soft X-ray flux and Ihe subsequent X-ray parameters have been computed from the 0.3-12 keV counts, where the source is detected. 
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The rest-frame 5 GHz flux density was computed from the 
FIRST (Becker, White, & Helfand 1995) or NVSS (Condon 
et al. 1998) observed 1.4 GHz flux density assuming a radio 
power-law slope of a = —0.8. The upper limits reported in 
the table are at the 3o level. The rest-frame 4400 A flux den- 
sity was computed from the AZ? 1450 ( 1+Z ) magnitude assuming 
an optical power-law slope of a = -0.79. Radio-loud quasars 
(RLQs) have R > 100, whereas radio-quiet quasars (RQQs) are 
characterized by R < 10 (e.g., Kellermann et al. 1989). 

SDSS 0913+5919 is the only radio-loud object (R w 41 1) in 
our sample and is also the highest redshift (z = 5.11) RLQ thus 
far detected in the X-ray band. For one quasar, BR 0305-4957, 
no radio data are available. None of the remaining quasars is a 
RLQ, although in some cases the observed-frame 1 .4 GHz flux- 
density constraints provided by the NVSS are not sufficiently 
tight to prove that a quasar is radio quiet (R < 10); at most 
these objects could be radio moderate (see Table 3). Given that 
w 85-90% of the quasar population is radio quiet (e.g., Stern 
et al. 2000), we expect the majority of the quasars with radio- 
loudness upper limits of 10-100 to be radio quiet. 

4. DISCUSSION 

Thanks mostly to Chandra (e.g., VOla; Brandt et al. 2002; 
Bechtold et al. 2003; V03) and ROSAT (e.g., KBS00; V03) 
observations, we have seen a substantial growth in the data 
needed to study the X-ray properties of z > 4 quasars. At 
present, there are « 70 quasars at z > 4 with X-ray detections. 10 
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FIG. 4. — Observed-frame, Galactic absorption-corrected 0.5-2 keV flux 
versus ABi 450 (i +z ) magnitude for z, > 4 AGNs. Object types are shown in the 
key. In particular, the objects presented in this paper are plotted as filled trian- 
gles (Chandra detections) and a filled square (XMM-Newton detection), while 
large, thick, downward-pointing arrows indicate X-ray upper limits for the ob- 
jects of our sample. Chandra RLQs are plotted as circled triangles. ROSAT 
upper limits (small downward-pointing arrows) are at the 3o confidence level. 
The slanted dashed lines show the a ox =— 1 .5 and a x=— 1-8 loci at z = 5.0 (the 
average redshift of the present sample). The dotted line shows the best-fit cor- 
relation reported in §4 [Equation (2)] for the whole sample of z, > 4 optically 
selected RQQs. 

Figure 4 shows the observed-frame, Galactic absorption- 
corrected 0.5-2 keV flux versus AB 1450 ( 1+ ,) magnitude for a 

10 See http://www.astro.psu.edu/users/niel/papers/highz-xray-detected.dat for a re: 



compilation of z > 4 Active Galactic Nuclei (AGNs). The de- 
tected objects presented in this paper are shown as filled trian- 
gles (10 Chandra detections) and a filled square (one XMM- 
Newton detection), while large and thick downward-pointing 
arrows indicate the two quasars in our sample with X-ray up- 
per limits. Because this figure is constructed from directly ob- 
served physical quantities, it is robust and useful for planning 
future X-ray observations of high-redshift quasars. We present 
a complementary discussion based on a ox , a parameter some- 
what less directly tied to observation, in §4.2. 

Quasars are well-known X-ray emitters at z ~ 0-3 (e.g., Vi- 
gnali et al. 1999; Reeves & Turner 2000; Page et al. 2003). 
Figure 4 indicates that X-ray emission is also a universal prop- 
erty of z > 4 quasars. 

Previous studies of z > 4 quasars have shown that there is 
a significant correlation between AZ? 1450 ( 1+Z ) and 0.5-2 keV 
flux (VOla; V03). The new X-ray detections and upper limits 
provided by the present work significantly improve the study 
of this correlation at fainter optical magnitudes. The average 
ABi45o(i +z ) magnitude of the sample presented here is 19.7, 
whereas the average for the bright PSS sample of V03 is 17.9. 
The presence of this correlation indicates that the same engine 
(namely accretion onto a supermassive black hole) is powering 
both the ultraviolet and X-ray emission. To study this corre- 
lation using all the available X-ray detections (41) and upper 
limits (52) for z > 4 optically selected RQQs (the majority of 
the objects plotted in Fig. 4), we have used the ASURV soft- 
ware package Rev 1.2 (La Valley, Isobe, & Feigelson 1992). 
To evaluate the significance of the correlation, we used sev- 
eral methods available in ASURV, namely the Cox regression 
proportional hazard (Cox 1972), the generalized Kendall's x 
(Brown, Hollander, & Korwar 1974), and the Spearman's p 
models. We have chosen to use only optically selected RQQs 
in order to have a homogeneous sample (see §4 of Vignali et 
al. 2002, hereafter V02, for the results of a sample of X-ray 
selected quasars; also see Silverman et al. 2002 and Castander 
et al. 2003). Once the Broad Absorption Line quasars (BALQ- 
SOs), which are known to be absorbed in the X-ray band (e.g., 
Brandt, Laor, & Wills 2000; Gallagher et al. 2001, 2002; Green 
et al. 2001; VOla), are excluded from the analysis (there are 
four in the original sample), we find that the AB^q^^ versus 
Fx correlation is significant at the > 99.9% level (all the tests 
above provide similar results). According to the EM (Estimate 
and Maximize) regression algorithm (Dempster, Laird, & Ru- 
bin 1977), the correlation is parameterized by 

log^x/erg cm- 2 s" 1 ) = [-(0.33 ±0.05) AS 



1450(l+z)" 



(7.89±0.91)], 
(2) 

which is plotted as a dotted line in Figure 4. Since 
ABi45o(i +z )=— 2.5 log / 1450 ^—48.6, the previous equation may 
also be written as 



-2 „-l\ 



log (F x /erg cm s 



(0.83 ±0.13) log/ 1450 . 



(3) 



In Figure 5 the quantity fx,obs/^x,fit is shown as a func- 
tion of AB 1450 (i +z ) magnitude for the z > 4 optically se- 
lected RQQs observed by Chandra. fx.obs represents the 
measured, Galactic absorption-corrected 0.5-2 keV flux, and 
Fxfit the flux expected from the best-fit correlation above. 
The majority of the z > 4 quasars observed by Chandra are 
within a factor of two of the best-fit correlation. The only 
quasar clearly showing a large difference between the mea- 
sured and expected fluxes is SDSS 0231—0728 (the bright- 

;ularly updated compilation of X-ray detections at z > 4. 
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est X-ray source in our sample). Since its X-ray emis- 
sion is unlikely to be contaminated by nearby X-ray sources 
(see Fig. 2), it appears plausible that its unusually large 
X-ray flux is due to an episode of variability. Unfortu- 
nately, the lack of an HET observation for SDSS 0231-0728 
does not allow us to examine this hypothesis further (as- 
suming that the X-ray and optical variability are correlated). 
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FIG. 5. — Plot of .Fx.obs/^x.fit as a function of ABi 450 ( 1+z ) magnitude for 
the z > 4 optically selected RQQs observed by Chandra. F x , bs represents the 
Galactic absorption-corrected 0.5-2 keV flux observed by Chandra and F x ,flt 
the flux expected on the basis of the best-fit correlation between the soft X-ray 
flux and AB 145 o( 1+z ) magnitude discussed in §4 [Equation (2)]. Symbols are 
the same as in Figure 4. The dotted line indicates where the measured and 
expected fluxes are the same. To avoid symbol crowding, we have only shown 
a representative average uncertainty on the flux ratios (the error bar near the 
middle right). 

A constant fit to the data points of Figure 5 (including the 
associated errors, but excluding the upper limits) is statistically 
unacceptable (% 2 = 54.6 for 33 degrees of freedom). Excluding 
SDSS 0231-0728 improves the quality of the fit (% 2 = 45.3), 
although it remains unacceptable. This indicates that there is 
a real dispersion in the quantity plotted in Figure 5 above the 
statistical noise (see also §4.2). 

4.1. X-ray Non-Detections: Likely Explanations 

All but two (SDSS 0338+0021 and SDSS 1737+5828) of 
the quasars presented here have been detected in the X-ray 
band. The non-detection of SDSS 0338+0021 is probably 
due to the short exposure time of the XMM-Newton observa- 
tion once the flaring-background intervals have been removed. 
Assuming the best-fit parameters of the AB 1450 (i +z ) versus 
0.5-2 keV flux correlation reported in Equation (2), we would 
expect a 0.5-2 keV flux of « 3.2 x 10~ 15 erg cm" 2 s" 1 . Al- 
though SDSS 0338+0021 was observed at the aimpoint, the on- 
axis flux limit of the XMM-Newton EPIC pn in a w 5.5 ks ex- 
posure (see Table 1) is a factor of ~ 2 higher than that achieved 
with Chandra in the same amount of time. 

No Chandra source counts were detected at the optical po- 
sition of SDSS 1737+5828. The faint ROSAT source reported 

11 This average value is flatter by 0.04 than that originally reported in V03 because 



as a possible counterpart of SDSS 1737+5828 by Anderson et 
al. (2001) is clearly present in the Chandra observation, but its 
X-ray position is 45" from the quasar's optical position. There 
are two possible explanations for the SDSS 1737+5828 non- 
detection: this quasar is either strongly absorbed or intrinsi- 
cally X-ray weak. If the non-detection is due to the presence 
of intrinsic X-ray absorption, a column density of Nh ^ 1 .2 x 
10 23 cm -2 is necessary to reproduce the Chandra constraints 
(see Tables 2 and 3) assuming the average a ox =— 1-71 of the 
sample (see §4.2) and an X-ray photon index of F = 2.0. The 
absorption possibility finds partial support from the published 
SDSS 1737+5828 optical spectrum (Anderson et al. 2001; see 
the grey curve in Fig. 3), which shows evidence of a strong, 
narrow N V absorption feature at z ~ 4.92. This absorption 
feature has a fairly high formal significance in the SDSS spec- 
trum, but given the importance of the intrinsic absorption inter- 
pretation, we wished to acquire an independent confirmation of 
the line; the new observation would also allow an investigation 
of any line variability. We therefore obtained a 30-min spec- 
troscopic observation of SDSS 1737+5828 with the HET (see 
§2.3). The resulting spectrum (black curve in Fig. 3) confirms 
the presence of a clear narrow N V absorption feature. This 
line is unresolved (the spectral resolution is m 1000 km s -1 ) 
and has a blueshift of « 700 km s -1 relative to the quasar red- 
shift of 4.94 reported by Anderson et al. (2001). In both the 
SDSS and HET spectra we see evidence for narrow Lya absorp- 
tion spanning the redshift of the N V absorption feature with 
the strongest absorption appearing at z ~ 4.90 and z ~ 4.945. 
There may be features at similar redshifts in Si IV and C IV, 
but the signal-to-noise ratio of the spectra prevents us from un- 
ambiguously identifying these systems. Nevertheless, the N V 
line alone strongly suggests that significant internal absorption 
is present in SDSS 1737+5828, since N V absorption of this 
strength is rarely caused by intervening galaxies (e.g., York et 
al. 1991) and the absorption is blueshifted with respect to the 
emission-line redshift. The Lya emission line appears weak; 
the rest-frame equivalent width is only m 35 A. On average, 
the rest-frame equivalent width of Lya + N V in high-redshift 
quasars is « 70-80 A (Schneider, Schmidt, & Gunn 1991; Fan 
et al. 2001). The evidence from the HET and SDSS spectra 
strongly suggests that significant internal absorption is present 
in SDSS 1737+5828, but the quasar is clearly not a strong high- 
redshift BALQSO such as SDSS 0856+5252 (z = 4.79; Ander- 
son et al. 2001), SDSS 1044-0125 (z = 5.74; Fan et al. 2000; 
Brandt et al. 2001; Maiolino et al. 2001; Goodrich et al. 2001), 
or SDSS 1 129-0142 (z = 4.85; Zheng et al. 2000). 

4.2. a ox Results 

Previous reports (e.g., VOla; Bechtold et al. 2003; V03) 
of steeper a ox indices for z ^ 4 quasars than for lower- 
luminosity, lower-redshift quasars are confirmed by the present 
analysis; see Figure 6. The 12 RQQs of our sample have 
(a ox }=— 1.71+0.05 (the error represents the standard deviation 
of the mean). This average a ox value is considerably steeper 
than that obtained for, e.g., the Bright Quasar Survey (BQS; 
Schmidt & Green 1983) RQQs atz < 0.5 ((a ox )=-1.56±0.02; 
see Brandt et al. 2000 and VOla). The sample of eight optically 
luminous z > 4 PSS RQQs presented in V03 has an even steeper 
average a ox ((a ox )=-1.77±0.03). u The PSS sample of V03 
has (log (Z2500 k/ er S s ~' Hz -1 )) = 32.2, while the one pre- 
sented here has (log (^500 A/ er 8 s ~' Hz -1 )} = 31.5. Accord- 

the flux correction discussed in §3. 
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ing to the known anti-correlation between a ox and 2500 A lu- 
minosity (e.g., Avni, Worrall, & Morgan 1995; Vignali, Brandt, 
& Schneider 2003, hereafter VBS03; also see the discussion 
later in this subsection) following the equation 12 

a ox = -(0.095 ±0.021) x / 250oA + (1.295 ±0.648), (4) 
we expect an average a ox of — 1 .70 for the sample of RQQs 
presented here, in good agreement with the average a ox value 
reported above. A summary of these results is shown in Fig- 
ure 6: the 0C ox distribution for the BQS RQQs (top panel) is 
clearly different from that of the z > 4 RQQs observed by 
Chandra thus far (bottom panel). In both oc ox distributions, 
some of the quasars with the steepest a ox values are BALQSOs 
or mini-BALQSOs; areas with arrows indicate upper limits. 



Sample Name 



BQS RQQk 



-±\ — i — I — i — h 



Chandra RQQs at zM 



N 

^detections 
dipper limits 


Bqs 
46 
45 
1 


16 
13 

3 


PSS+BK1 
19 
18 
1 


Range of z 
Median z. 


0.061-0.472 
0.17±0.01 
0.16 


4.62-6.28 
5.11±0.11 
4.96 


4.00^.73 
4.28±0.04 
4.23 


Range of log / 
log I 

* & 2500 A ' 

Median los / 

6 2500 A 


29.14-30.70 


31.22-31.75 


31.54-32.59 


29.66±0.05 


31.53±0.04 


31.96±0.06 


29.60 


31.55 


31.92 


Range of log L 2 -iokeV 

{ log Z-2-lOkcV > " 

Median log L 2 _ 10 kcV 


43.36^14.79 
44.07±0.10 
44.13 


44.53^15.89 
44.99±0.09 
45.03 


44.75^15.82 
45.34±0.07 
45.41 


Range of a ox 
{ a,« > * 
Median a ox 


-1.23 to -2.06" 
-1.56 ±0.02 
-1.55 


-1.43 to -1.82 
-1.69 ±0.03 
-1.71 


-1.56 to -1.94 
-1.73 ±0.02 
-1.73 


References 


(1,2,3,4) 


(4,5,6,7,8) 


(5,7,9) 



NOTE. — The luminosities are in the quasars' rest frames and have been 



obtained (or converted) adopting Hq=70 km s Mpc 



in a A-cosmology 



-2.5 



1.5 



with £2m=0.3 and £2a=0.7. / ^ is in units of erg s 1 Hz 1 and L2-10 keV 

is in units of erg s~ 1 . Corrections for the Chandra ACIS quantum efficiency 
decay at low energies have been applied in the computation of the X-ray 
luminosities and a ox values (see §3 for details). 

BALQSOs and mini-BALQSOs have been excluded from all of the sam- 
ples used here. For the SDSS and PSS+BRI samples at z > 4, the values 
presented here have been obtained using Chandra and XMM-Newton re- 
sults. For the BQS sample, the rest-frame 2-10 keV luminosities have been 
obtained from a representative subsample of 21 BQS RQQs observed by 
ASCA (George et al. 2000) and BeppoSAX (Mineo et al. 2000) in the ob- 
served ks 0.5-10 keV and ss 0.1-10 keV bands, respectively. This subsam- 
ple covers approximately the same range in redshift as the original sample 
of 46 BQS RQQs with M B < -23 used in §3.1.1 of Brandt et al. (2001). 
a The errors on average values represent the standard deviation of the mean. 
b ct ox =-2.06 represents the slope of PG 2214+139. See Appendix A of 
Laor & Brandt (2002) for a discussion about possible absorption in this 
quasar. 

REFERENCES. — (1) Brandt et al. 2000; (2) George et al. 2000; (3) Mineo 
et al. 2000; (4) Brandt et al. 2001; (5) VOla; (6) Brandt et al. 2002; (7) 
Bechtold et al. 2003; (8) this paper; (9) V03. 



FIG. 6.— a ox distributions for the BQS RQQs observed by ROSAT (top) 
and z > 4 RQQs observed by Chandra (bottom). The shaded areas indicate 
the BALQSOs or mini-BALQSOs; the areas with arrows indicate 95% confi- 
dence upper limits. The average uncertainty in the a ox measurements from the 
Chandra observations is shown in the bottom panel as the error bar near the 
upper right. 

The principal properties (redshifts, rest-frame 2500 A and 
2-10 keV luminosities, and a ox values) of representative sam- 
ples of RQQs at low and high redshifts are reported in Table 4. 

We find no significant correlation of a ox with either rest- 
frame 2500 A luminosity or redshift using only the available 
z > 4 RQQs. However, we can also use the new X-ray detec- 
tions and upper limits presented here to extend the VBS03 study 
of a ox dependence upon 2500 A luminosity and redshift. 13 As 
in VBS03, we mostly use the SDSS RQQs taken from the Early 
Data Release (EDR) quasar catalog (Schneider et al. 2002). In 
particular, at z < 4 we use the 128 EDR RQQs of VBS03 (see 
their §3.1.1) covered by pointed ROSAT PSPC and HRI obser- 
vations, with an overall X-ray detection fraction of « 48%. At 
z > 4 we use all published Chandra (VOla; Brandt et al. 2002; 
Bechtold et al. 2003) and XMM-Newton (Brandt et al. 2001) 
observations of SDSS RQQs in addition to the ones presented 
in this paper. 

TABLE 4 

Principal Properties of Three Optically Selected Samples 
of rqqs at low and high redshifts 

12 This is an updated version of Equation (2) of VBS03, taking into account the correction for the quantum efficiency decay of ACIS. 

13 See §2.2 of VBS03 for discussion about the method used to derive 2500 A fluxes from the SDSS filters. 



In total we have 20 SDSS RQQ observations at z > 4 (includ- 
ing RD 657, observed in one of the fields surveyed by the 
SDSS; see §2.2 and Djorgovski et al. 2003), 16 of which 
have X-ray detections. To improve the statistics at high lu- 
minosities and redshifts without introducing any evident bias 
(see §3.2 of VBS03), we also include 15 PSS RQQs and seven 
BRI RQQs at z > 4 observed by Chandra (VOla; Bechtold 
et al. 2003; V03; Priddey et al., in preparation); only one 
of these is not detected in the X-ray band (VOla). There- 
fore, the final sample used in the present analysis is comprised 
of 170 RQQs (w 58% with X-ray detections) over the red- 
shift range z ~ 0.2-6.3 and the rest-frame 2500 A luminosity 
range log (/2500 k/ er S s 1 Hz -1 ) « 29.0-32.6. Eleven of these 
quasars are known BALQSOs; results obtained both including 
and excluding BALQSOs are presented below. As luminosity 
is correlated with redshift in flux-limited samples, we use a par- 
tial correlation analysis technique developed for use with cen- 
sored data (Akritas & Siebert 1996). Partial correlation analy- 
ses confirm previous findings (e.g., Avni et al. 1995; VBS03): 
0C ox depends primarily upon rest-frame 2500 A luminosity. The 
significance of this anti-correlation is 3.5^-.0o (slightly higher 
than in VBS03), depending upon whether PSS/BRI quasars and 
BALQSOs are included in the analyses (see Table 5). 



TABLE 5 

DEPENDENCIES OF (X„ x UPON REDSHIFT AND 2500 A LUMINOSITY DENSITY 
FROM PARTIAL CORRELATION ANALYSIS 
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Sample Name 


Number 


Oiox-Z 


° x 6 2500 




of RQQs 




SDSS only, with BALQSOs 


148 


1.2a 


3.7a 


SDSS only, without BALQSOs 


138 


0.9a 


3.5a 


SDSS + PSS, with BALQSOs 


163 


1.3a 


3.9a 


SDSS + PSS, without BALQSOs 


152 


1.2a 


3.9a 


SDSS + PSS + BRI, with BALQSOs 


170 


1.4a 


4.0a 


SDSS + PSS + BRI, without BALQSOs 


159 


1.1a 


4.0a 



No significant correlation of a ox with redshift is found. These 
results suggest that the mechanism driving quasar broad-band 
emission is similar in the local and early Universe, with no 
evidence for unusual phenomena in high-redshift quasar nu- 
clei such as "trapping radius" effects (e.g., Begelman 1979) or 
accretion-disk instabilities (e.g., Lightman & Eardley 1974). 

Although we did not find a significant correlation of a ox with 
either rest-frame 2500 A luminosity or redshift using only the 
Z > 4 RQQs, the spread of a ox for z > 4 RQQs is unlikely to be 
explained solely by measurement errors. To demonstrate this, 
we have applied the likelihood method described in Section IV 
of Maccacaro et al. (1988) to calculate the intrinsic dispersion 
of a ox (accounting for measurement errors). Ignoring upper 
limits, we find an intrinsic dispersion of O; = O.O62+QQ33 (68% 
confidence errors for two interesting parameters; see Maccac- 
aro et al. 1988). Accounting for upper limits would only in- 
crease the intrinsic dispersion (see the bottom panel of Fig. 6). 
For comparison, the "straight" dispersion (not accounting for 
measurement errors) is o s = 0.1 10. 

SDSS 023 1 -0728 and RD 657 have the flattest a ox values in 
our sample. It is possible (see §4) that they experienced vari- 
ability over the period between the optical and X-ray observa- 
tions. An alternative explanation could be the presence of some 
dust obscuring the optical/ultraviolet radiation but not the X-ray 
emission. The optical spectra of these quasars, however, do not 
show any clear evidence supporting this scenario. Radio emis- 
sion can affect the quasar broad-band spectral energy distribu- 
tions producing flatter a ox values. Although significant radio 
emission can be ruled out for SDSS 0231—0728 (see Table 3), 
the radio-loudness upper limit for RD 657 is not sufficiently 
tight to exclude some radio emission from this quasar. 

SDSS 0913+5919, the only radio-loud object in our sam- 
ple and the highest redshift RLQ with an X-ray detection, 
has the highest R parameter (R « 411) and the steepest a ox 
( — 1 .67 ± 0. 1 1) among the RLQs observed thus far by Chandra 
atz >4 (VOla; V03). X-ray absorption might be an explanation 
for the steep a ox value, although the detection of five of the six 
source counts in the soft band (see Table 2) suggests a different 
interpretation. The anti-correlation between a ox and rest-frame 
2500 A luminosity (e.g., Avni et al. 1995; VBS03) cannot ex- 
plain its steep (for a RLQ) a ox value, since SDSS 0913+5919 
is also the optically faintest z > 4 RLQ observed by Chandra to 
date. 

4.3. Joint X-ray Spectral Fitting of SDSS Quasars at z> 4.8 

To define the average X-ray spectral properties of high- 
redshift SDSS quasars and to make a comparison with those 
obtained for a sample of nine optically luminous PSS quasars 
at z« 4.1-4.5 (V03), we have selected a sample of 13 SDSS 
quasars at z > 4.8 detected by Chandra (( z ) = 5.23). Their 
total exposure time is 81.6 ks. The quasars in this sample are 
selected to have more than two counts in the full band; this 

If we increase the number of soft-band counts by 10% to take into account crudely the quantum efficiency decay of ACIS at soft energies, we obtain r=1.93±g;^. 
The small differences in the photon indices obtained via joint spectral fitting and band-ratio analysis (with a 10% increment in the number of soft-band counts) are 
probably due to the fact that only the former retains all the spectral information. 



choice excludes SDSS 0756+4502 (presented in this paper) 
and SDSS 1208+0010 (VOla). Five of the quasars used in 
this sample have been previously published: the three highest 
redshift SDSS quasars observed thus far (SDSS 1030+0524, 
SDSS 1306+0356, and SDSS 0836+0054; Brandt et al. 
2002); SDSS 1204-0021 (Bechtold et al. 2003); and 
SDSS 0211-0009 (VOla). Eleven of the 13 quasars are ra- 
dio quiet; the exceptions are SDSS 0836+0054 (radio moderate 
with R « 10) and SDSS 0913+5919 (radio loud with R w 41 1; 
see Table 3). In the following we present the joint X-ray spec- 
tral analysis both with and without RLQs, since RLQs often 
appear to be characterized by flatter X-ray photon indices (e.g., 
Cappi et al. 1997). The sample does not appear to be bi- 
ased by the presence of a few high signal-to-noise ratio ob- 
jects. For example, the removal of the radio-moderate quasar 
SDSS 0836+0054, the second brightest source in the present 
sample, does not provide significantly different results (see be- 
low). 

Source counts were extracted from 2" radius circular aper- 
tures centered on the X-ray position of each quasar. The back- 
ground was taken from annuli centered on the quasars, avoid- 
ing the presence of nearby X-ray sources. In total we have 120 
counts (in the rest-frame w 2.5-36 keV band) from the source 
extraction regions and five counts from the background regions 
normalized to the size of those used for the sources; half of the 
background counts are from the SDSS 0941+5947 field (see 
§2.1). All of the ARFs used in the spectral analysis have been 
corrected to account for the low-energy quantum efficiency de- 
cay of ACIS (see §2.1). As in §2.1, spectral analysis was car- 
ried out with XSPEC using the background-subtracted data and 
C-statistic (Cash 1979). Since C-statistic fitting is performed 
using the unbinned data, it retains all the spectral information 
and allows one to associate with each quasar its own Galactic 
absorption column density and redshift (for the fitting with in- 
trinsic absorption). 

Joint spectral fitting with a power-law model (leaving the 
normalizations free to vary) and Galactic absorption provides 
a good fit to the data as judged by the small data-to-model 
residuals and checked with 10,000 Monte-Carlo simulations. 
The resulting photon index is Y = 1.84^1^. Excluding the 
two non-RQQs from the joint spectral fitting gives a similar 
result (r = 1.86^03'). These photon indices are consistent 
with those obtained for z w 0-3 samples of RQQs in the rest- 
frame 2-10 keV band (e.g., George et al. 2000; Mineo et 
al. 2000; Reeves & Turner 2000; Page et al. 2003) and the 
z « 4.1^.5 optically luminous PSS quasars presented in V03 
(r = 1.98 ±0.16). As a consistency check, we compared the 
photon index obtained from joint spectral fitting with that de- 
rived from band-ratio analysis (i.e., the ratio of the 2-8 keV to 
0.5-2 keV counts; we have assumed for the Galactic column 
density the average value for the 13 quasars weighted by their 
number of counts). In this case we obtain r = 1 .85^ 20- 14 

We also constrained the presence of neutral intrinsic ab- 
sorption; the joint-fitting technique described above provides 
a counts-weighted average column density. Solar abundances 
have been assumed, although previous optical studies indicate 
that high-redshift quasar nuclei are often characterized by su- 
persolar abundances of heavy elements (e.g., Hamann & Fer- 
land 1999; Constantin et al. 2002; Dietrich et al. 2003). Dou- 
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bling the abundances in the fit gives a reduction in the column 
density of a factor of w 2. No evidence for absorption in ex- 
cess of the Galactic column density has been found, as shown 
by the 68, 90, and 99% confidence regions plotted in Figure 7. 
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FIG. 7. — Confidence regions for the photon index and intrinsic column den- 
sity derived from joint spectral fitting of the 13 SDSS quasars (11 RQQs, two 
non-RQQs) at z > 4.8 thus far detected by Chandra with more than two full- 
band counts (solid contours; see §4.3 for details). Dashed contours indicate 
the results obtained when only the 11 SDSS RQQs at z > 4.8 are considered 
in the joint spectral fitting. 

On average, it appears that any neutral X-ray absorption 
in z > 4.8 X-ray detected SDSS quasars has a column den- 
sity of N H < 3.95 x 10 22 cm" 2 (at the 90% confidence level; 
N H < 4.37 x 10 22 cirr 2 when only the RQQs are used in 
the spectral analysis). For comparison, X-ray absorption of 
ss (2-5) x 10 22 cm~ 2 has been detected in a few RQQs and 
some RLQs at z ~ 1-3 (e.g., Cappi et al. 1997; Reeves & 
Turner 2000). Given the limited counting statistics available 
for the present sample, it is not possible to obtain useful spectral 
constraints for more complex absorption models (e.g., ionized 
absorption). The typical upper limits on the equivalent widths 
of narrow iron Ka lines (either neutral or ionized) are in the 
range 420-950 eV. 

The X-ray spectral results at z > 4 presented here and in 
V03 allow us to place constraints on any spectral changes in 
the X-ray continua of RQQs as a function of redshift. In Fig- 
ure 8 we plot the photon index as a function of redshift for some 
optically selected samples of RQQs. The photon indices plot- 
ted in the figure are for the nuclear, high-energy ( ^ 2 keV) 
component, in order to avoid spectral complexities such as 
soft excesses. The ASCA RQQs have been taken from Vig- 
nali et al. (1999; V99 in the key), George et al. (2000; GOO), 
Reeves & Turner (2000; RT00), and Vignali et al. (2001b; 
VOlb). The BeppoSAX RQQs are from Mineo et al. (2000; 
MOO). At the highest redshifts we plot the z = 3.91 gravitation- 
ally lensed BALQSO APM 08279+5255 (Chartas et al. 2002; 
C02), the results obtained from joint X-ray spectral fitting of 
nine PSS quasars in the redshift range z ~ 4.1-4.5 (V03), and 
the 1 1 SDSS RQQs at z > 4.8 presented in this paper. The 
distribution of X-ray photon indices is characterized by signif- 
icant scatter, so precise parametric modeling of any F versus 
redshift relation is difficult. However, according to a Spear- 
man's p correlation test, a correlation between F and redshift 



is not significant (2.3a). Applying a linear regression algo- 
rithm that accounts for intrinsic scatter and measurement errors 
(Akritas & Bershady 1996), we obtain | |f | < 0.05 ±0.03. 
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FIG. 8. — Plot of photon index versus redshift for optically selected RQQs. 
The ASCA RQQs have been taken from Vignali et al. (1999; V99 in the key), 
George et al. (2000; GOO), Reeves & Turner (2000; RT00), and Vignali et al. 
(2001b; VOlb). The BeppoSAX RQQs are from Mineo et al. (2000; MOO). 
At the highest redshifts we plot the z = 3.91 gravitationally lensed BALQSO 
APM 08279+5255 (Chartas et al. 2002; C02), the results obtained from joint 
X-ray spectral fitting of nine PSS quasars in the redshift range z ~ 4.1-4.5 
(V03), and the 11 SDSS RQQs at z > 4.8 presented in this paper (the latter 
two data points are plotted at the average redshifts of the relevant papers). Er- 
ror bars are at the 68% confidence level for all data points. The object with the 
steepest X-ray photon index in this plot, PG 1543+489 (George et al. 2000), 
is known to be a narrow-line type 1 quasar (Boroson & Green 1992). 



5. SUMMARY 

We have reported the Chandra and XMM-Newton observa- 
tions of a sample of 13 quasars at z ~ 4.7-5.4, the majority of 
them selected by the SDSS, and compared their properties with 
those of all the z > 4 quasars previously observed in the X-ray 
band. The principal results are the following: 

• Eleven quasars have been detected in the X-ray band, 
doubling the number of z > 4.8 X-ray detected quasars. 
The X-ray detections include SDSS 0913+5919, the 
highest redshift (z = 5.11) RLQ thus far detected in 
the X-ray band. Two quasars, SDSS 0338+0021 and 
SDSS 1737+5828, were not detected by the present 
observations. The former non-detection is probably 
due to the short exposure time of the XMM-Newton 
observation after removal of flaring-background inter- 
vals. For the latter quasar, strong X-ray absorption 
(7Y H ^l-2xl0 23 cm- 2 at the source redshift) is a likely 
explanation. An HET spectrum of this source, revealing 
absorption features, seems to support this hypothesis. 
For our sample, no evidence for source extension or an 
overdensity of companion objects has been found. 

• The presence of a significant correlation between the 
soft X-ray flux and AB 1450 (i +z ) magnitude has been 
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extended to the highest redshifts. This indicates that 
the same engine (namely accretion onto a supermassive 
black hole) is powering the ultraviolet and X-ray emis- 
sion. 

• The anti-correlation between a ox and rest-frame 2500 A 
luminosity density found for SDSS RQQs over the red- 
shift range z ~ 0.2-6.3 (VBS03) has been confirmed 
(3.5-4.0O significance level) after the inclusion of the 
RQQs presented in this paper. Similarly to VBS03, 
no significant correlation between a ox and redshift has 
been found. These results suggest that the mechanism 
driving quasar broad-band emission is similar in the lo- 
cal and early Universe, with no evidence for unusual 
phenomena in high-redshift quasar nuclei such as "trap- 
ping radius" effects or accretion-disk instabilities. 

• The joint w 2.5-36 keV rest-frame X-ray spectrum of 
the 13 SDSS quasars at z > 4.8 observed thus far by 
Chandra is well parameterized by a simple power-law 
model with Y = 1.84+Jj ^ and no evidence for intrinsic 
absorption (Nn £ 4.0 x 10 22 cirr 2 at 90% confidence). 

• There is no evidence for significant (> 30%) optical 
variability over the time interval of a few years (in the 
observed frame) between the SDSS and HET observa- 
tions. 

Among the quasars presented in this paper, the only 
object suitable for moderate-quality (m 1000 photons) 
X-ray spectroscopy in a 50 ks XMM-Newton observation 
is SDSS 0231-0728. X-ray spectroscopy of the other 
quasars will be possible only with future missions such as 
Constellation-X, XEUS, and Generation-X (see V03 for some 
relevant simulations). 
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APPENDIX 

Here we report the updated X-ray flux values (and derived quantities) for z > 4 quasars observed by Chandra in the past. The 
X-ray values have been corrected for the ACIS quantum energy decay at low energies. 



Table A 1 

Optical, X-ray, and Radio Properties of z > 4 Quasars Previously Observed by Chandra 



Object 


z 


AB 1450(l+z) 


■^2500 A 




M B 


A.5-2 kcV 


fl keV 


log(L 2 kcv) b 


lOg(L 2 _10kev) 




R 


PSS 0059+0003 


4.16 


19.3 


1.06 


31.54 


-27.62 


10.2 


7.88 


27.41 


45.30 


-1.59 


<3.3 


BR 0103+0032 


4.43 


18.7 


1.85 


31.82 


-28.33 


20.1 


16.28 


27.77 


45.66 


-1.56 


2.9 


PSS 0121+0347 


4.13 


18.5 


2.22 


31.86 


-28.41 


41.3 


31.65 


28.01 


45.90 


-1.48 


300.3 


PSS 0133+0400 


4.15 


18.2 


2.93 


31.98 


-28.72 


20.8 


16.01 


27.72 


45.61 


-1.64 


<4.4 


PSS 0134+3307 


4.53 


18.5 


2.22 


31.92 


-28.57 


11.5 


9.45 


27.55 


45.44 


-1.68 


<6.1 


PSS 0209+0517 


4.14 


17.8 


4.24 


32.14 


-29.12 


15.2 


11.66 


27.58 


45.47 


-1.75 


<3.0 


SDSS 0210-0018 


4.77 


19.3 


1.06 


31.63 


-27.85 


37.3 


32.12 


28.11 


46.00 


-1.35 


86.1 


SDSS 0211-0009 


4.90 


20.0 


0.56 


31.37 


-27.20 


3.8 


3.36 


27.15 


45.04 


-1.62 


<7.5 


BR 0241-0146 


4.06 


18.4 


2.44 


31.88 


-28.48 


4.3 


3.26 


27.01 


44.84 


-1.87 


< 1.5 


PSS 0248+1802 


4.43 


18.1 


3.21 


32.06 


-28.93 


21.2 


17.16 


27.79 


45.68 


-1.64 


<4.2 


BR 0308-1734 


4.00 


18.1 


3.21 


31.99 


-28.76 


4.1 


3.06 


26.97 


44.86 


-1.93 


< 3.9 


BR 0401-1711 


4.23 


18.9 


1.54 


31.71 


-28.05 


12.9 


10.07 


27.53 


45.42 


-1.61 


< 8.5 


SDSS 0836+0054 


5.82 


18.8 


1.67 


31.96 


-28.67 


12.3 


12.50 


27.83 


45.73 


-1.58 


9.8 


PSS 0926+3055 


4.19 


16.7 


11.67 


32.59 


-30.24 


33.5 


25.94 


27.93 


45.82 


-1.79 


<0.3 


PSS 0955+5940 


4.34 


18.4 


2.44 


31.93 


-28.60 


6.0 


4.79 


27.22 


45.12 


-1.81 


<5.4 


PSS 0957+3308 


4.20 


18.2 


2.93 


31.99 


-28.74 


11.3 


8.78 


27.46 


45.36 


-1.74 


< 1.1 


SDSS 1030+0524 


6.28 


19.7 


0.76 


31.67 


-27.94 


2.5 


2.73 


27.22 


45.12 


-1.71 


<6.8 


BR 1033-0327 


4.51 


18.8 


1.69 


31.80 


-28.26 


8.4 


6.87 


27.41 


45.30 


-1.69 


<2.5 


PSS 1057+4555 


4.12 


17.6 


5.09 


32.22 


-29.31 


25.3 


19.34 


27.80 


45.68 


-1.70 


2.3 


SDSS 1129-0142 


4.85 


19.2 


1.17 


31.68 


-27.98 


<2.4 


< 1.75 


< 26.86 


< 44.84 


< -1.85 


<3.4 


SDSS 1204-0021 


5.03 


19.1 


1.28 


31.75 


-28.14 


4.6 


4.13 


27.26 


45.15 


-1.72 


< 3.4 


SDSS 1208+0010 


5.27 


20.5 


0.35 


31.22 


-26.82 


1.3 


1.24 


26.77 


44.66 


-1.71 


< 13.1 


PC 1247+3406 


4.90 


19.2 


1.17 


31.69 


-28.00 


11.0 


9.66 


27.61 


45.50 


-1.57 


< 3.3 


SDSS 1306+0356 


5.99 


19.6 


0.85 


31.68 


-27.97 


5.7 


5.94 


27.53 


45.42 


-1.59 


<5.3 


PSS 1317+3531 


4.36 


18.9 


1.54 


31.73 


-28.10 


3.6 


2.91 


27.01 


44.90 


-1.81 


<2.3 


PSS 1326+0743 


4.09 


17.2 


7.36 


32.37 


-29.70 


29.4 


22.32 


27.85 


45.74 


-1.73 


<0.5 


PSS 1347+4956 


4.51 


17.4 


6.12 


32.36 


-29.66 


17.5 


14.35 


27.73 


45.62 


-1.78 


0.1 


PSS 1435+3057 


4.35 


19.1 


1.28 


31.65 


-27.90 


<3.6 


<2.87 


< 27.00 


< 44.90 


< -1.79 


<2.9 


PSS 1443+2724 


4.42 


19.0 


1.40 


31.70 


-28.03 


11.3 


9.10 


27.51 


45.41 


-1.61 


<2.8 


PSS 1443+5856 


4.26 


17.7 


4.64 


32.20 


-29.26 


4.7 


3.66 


27.09 


44.99 


-1.96 


<2.8 


SDSS 1532-0039 


4.62 


19.4 


0.97 


31.57 


-27.70 


<2.5 


<2.10 


< 26.91 


< 44.80 


< -1.79 


<4.3 


SDSS 1605-0112 


4.92 


19.4 


0.97 


31.61 


-27.80 


<3.0 


<2.62 


< 27.05 


< 44.94 


< -1.75 


<4.5 


BR 2212-1626 


4.00 


18.6 


2.03 


31.79 


-28.25 


13.4 


9.95 


27.48 


45.38 


-1.65 


<6.2 



Note. — Units are the same as in Table 3. 
"Rest-frame 2500 A luminosity density (erg s _1 Hz -1 ). 
b Rest-frame 2 keV luminosity density (erg s _1 Hz -1 ). 
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Chandra AND XMM-Afevvfon OBSERVATIONS OF THE FIRST QUASARS: X-RAYS FROM THE AGE OF 

COSMIC ENLIGHTENMENT 

C. Vignali, ' W. N. Brandt, 1 D. P. Schneider, 1 S. F. Anderson, 2 X. Fan, 3 J. E. Gunn, 4 S. Kaspi, 5 

G. T. Richards, 4 and Michael A. Strauss 4 
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ABSTRACT 

We report on Chandra and XMM-Newton observations of a sample of 13 quasars at z « 4.7-5.4 mostly taken 
from the Sloan Digital Sky Survey (SDSS). The present sample complements previous X-ray studies of z > 4 
quasars, in which the majority of the objects are optically more luminous and at lower redshifts. All but two of 
our quasars have been detected in the X-ray band, thus doubling the number of z > 4.8 X-ray detected quasars. 
The two non-detections are likely to be due to a short exposure time (SDSSp J033829. 31+002156. 3) and to the 
presence of intrinsic absorption (SDSSp J173744. 87+582829. 5). We confirm and extend to the highest redshifts 
the presence of a correlation between AB 1450(1 +z ) magnitude and soft X-ray flux for z > 4 quasars, and the presence 
of a steeper optical-to-X-ray spectral energy distribution (parameterized by a ox ) for high-luminosity, high-redshift 
quasars than for lower-luminosity, lower-redshift quasars. The second effect is likely due to the known anti- 
correlation between oc ox and rest-frame 2500 A luminosity, whose significance is confirmed via partial correlation 
analysis. The joint w 2.5-36 keV rest-frame spectrum of the z > 4.8 SDSS quasars observed thus far by Chandra 
is well parameterized by a power-law with photon index T= 1 .84^q lh this photon index is consistent with those of 
z ~ 0-3 quasars and that obtained from joint spectral fitting of z ~ 4.1-4.5 optically luminous Palomar Digital Sky 
Survey quasars. No evidence for widespread intrinsic X-ray absorption has been found (Na ~ 4.0 x 10 22 cm~ 2 
on average at 90% confidence). We also obtained Hobby-Eberly Telescope (HET) photometric observations for 
eight of our target quasars. None of these shows significant (> 30%) optical variability over the time interval of a 
few years (in the observed frame) between the SDSS and HET observations. 



1. INTRODUCTION 

One of the main themes in modern astronomy is the study 
of the earliest massive objects to form in the Universe, now 
known as far back as the end of the reionization epoch (e.g., 
Rees 1999; Loeb & Barkana 2001; Fan et al. 2002; Hu et al. 
2002). Quasars at z 4 can be studied in detail; because of their 
large broad-band luminosities, they are typically detectable at 
most wavelengths [e.g., in the radio, Schmidt et al. 1995; Stern 
et al. 2000; in the millimeter, Omont et al. 2001; Carilli et al. 
2001; in the sub-millimeter, McMahon et al. 1999; Isaak et al. 
2002; in the near-infrared, Barkhouse & Hall 2001; Pentericci 
et al. 2003; in the optical, Schneider, Schmidt, & Gunn 1989; 
Fan et al. 2001; in the X-ray band, Kaspi, Brandt, & Schneider 
2000 (hereafter KBS00); Vignali et al. 2001, 2003 (hereafter 
VOla, V03); Bechtold et al. 2003]. Ground-based, large-area 
optical surveys such as the Sloan Digital Sky Survey (SDSS; 
York et al. 2000), the Automatic Plate Measuring facility sur- 
vey (BRI; Irwin, McMahon, & Hazard 1991), and the Palomar 
Digital Sky Survey (PSS; Djorgovski et al. 1998) have discov- 
ered more than 300 quasars at z > 4 thus far up to z = 6.43 (Fan 
et al. 2003). 6 

To extend our knowledge of the X-ray properties of z > 4 
quasars, we have started a program to observe, with Chan- 
dra and XMM-Newton, both the optically most luminous 
(M B « -28.4 to -30.2) z~ 4.1-4.5 PSS/BRI quasars (e.g., 

1 Department of Astronomy and Astrophysics, The Pennsylvania State University, 525 Davey Laboratory, University Park, PA 16802, USA; chris@astro.psu.edu, 
niel@astro.psu.edu, and dps@astro.psu.edu. 

2 Astronomy Department, University of Washington, Box 351580, Seattle, WA 98195-1580, USA; anderson@astro.washington.edu. 



V03) and the higher-redshift, optically fainter SDSS quasars 
(VOla; Brandt et al. 2002; see Brandt et al. 2003 for a recent 
review). Selecting z > 4 quasar samples for X-ray observations 
in two different regions of the luminosity-redshift parameter 
space allows study of any dependence of the X-ray properties 
upon the optical selection criteria. In particular, it is possible to 
test for dependencies of z > 4 quasar X-ray continua on optical 
luminosity or redshift. 

Here we report the results obtained from 1 1 Chandra obser- 
vations and one XMM-Newton observation of a sample of 13 
high-redshift (z « 4.7-5.4), moderately luminous (Mb ~ —26.1 
to —28.5) quasars, mostly discovered by the SDSS and pub- 
lished in Anderson et al. (2001). One of these quasars (RD 657; 
Djorgovski et al. 2003) was serendipitously observed and de- 
tected in the XMM-Newton observation (i.e., there were two 
objects in one field). The only target quasar not coming from 
the SDSS data is BR B0305-4957 (Storrie-Lombardi et al. 
2001; Peroux et al. 2001). Given that the accretion rate 
(relative to the Eddington rate) of quasars may increase with 
redshift (e.g., Kauffmann & Haehnelt 2000), the first quasars 
could show unusual phenomena such as "trapping radius" ef- 
fects (e.g., Begelman 1979) or accretion-disk instabilities (e.g., 
Lightman & Eardley 1974). Furthermore, our targets are bright 
enough to assure in most cases X-ray detections with snap- 
shot (e.g., « 4-10 ks) observations, given the correlation found 



3 Steward Observatory, The University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA; fan@as.arizona.edu. 

4 Princeton University Observatory, Peyton Hall, Princeton, NJ 08544-1001, USA; jeg@astro.princeton.edu, 



gtr@astro.princeton.edu, 
strauss@astro.princeton.edu. 

5 School of Physics and Astronomy, Raymond and Beverly Sackler Faculty of Exact Sciences, Tel-Aviv University, Tel-Aviv 69978, Israel; shai@wise.tau.ac.il. 
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and 



' See http://www.astro.caltech.edu/~george/z4.qsos for a listing of known high-redshift quasars. 
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between AB 1450 ( 1+?; ) magnitude and soft X-ray flux (see, e.g., 
VOla; V03). The number of z > 4.8 X-ray detected quasars is 
approximately doubled by the objects presented here. As shown 
in Figure 1, where the AB 1450 ( 1+Z \ magnitude is plotted as a 
function of redshift, the objects presented in this paper suitably 
complement previous X-ray studies of z > 4 quasars. All but 
one of these quasars are likely radio-quiet (see §3 for details). 




4 4.5 5 5.5 6 

Redshifl 

FIG. 1. — AB 1450 (i +z ) versus redshift for the known z > 4 AGNs [small 
dots; taken from Djorgovski as of 2002 July (see Footnote 6)]. The X-ray de- 
tected (undetected) quasars presented in this paper are shown as filled (open) 
triangles (Chandra observations) and a filled (open) square (XMM-Newton 
observation). Filled and open circles indicate z > 4 X-ray detected AGNs and 
tight upper limits, respectively, from previous observations [Schneider et al. 
1998; Waddington et al. 1999 (note that this identification is still debated; 
L. L. Cowie 2002, private communication); KBS00; VOla; Brandt et al. 2001, 
2002; Silverman et al. 2002; Vignali et al. 2002 (V02); Barger et al. 2002; 
Bechtold et al. 2003; V03; Castander et al. 2003], most of which were per- 
formed by Chandra and ROSAT. 

Hereafter we adopt Hq=70 km s -1 Mpc -1 in a A-cosmology 
with i2 M =0.3 and il A =0.7 (e.g., Turner 2002). 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Chandra Observations 

Eleven of the 13 quasars were observed by Chandra dur- 
ing Cycle 3 with « 3.7-9.8 ks observations (see Table 1; here- 
after all of the quasars will be referred to via their abbreviated 
names). All of the sources were observed with the Advanced 
CCD Imaging Spectrometer (ACIS; Garmire et al. 2003) with 
the S3 CCD at the aimpoint. Faint mode was used for the 
event telemetry format, and ASCA grade 0, 2, 3, 4 and 6 events 
were used in the analysis. This particular grade selection is 
standard and appears to optimize the signal-to-background ra- 
tio (see §6.3 of the Chandra Proposers' Observatory Guide). 
Source detection was carried out with WAVDETECT (Freeman 
et al. 2002). We used wavelet transforms (using a Mexican hat 
kernel) with wavelet scale sizes of 1, 1.4, 2, 2.8, and 4 pixels 
and a false-positive probability threshold of 10~ 4 . Given the 
small number of pixels being searched due to the known source 

7 To account for the quantum efficiency decay of ACIS at low energies, possibly 
dependent correction to the ACIS quantum efficiency. We generated new ancilk 
throughout this paper. 



positions and the excellent angular resolution of Chandra, the 
probability of spurious detections is extremely low; most of 
the sources were in fact detected at a false-positive probabil- 
ity threshold of 10~ 6 . Typically, detections were achieved with 
a wavelet scale of 1.4 pixels or less. 

Source searching was performed in four energy ranges for 
consistency with our previous work (VOla; V03): the ultra- 
soft band (0.3-0.5 keV), the soft band (0.5-2 keV), the hard 
band (2-8 keV), and the full band (0.5-8 keV); in the redshift 
range of z ~ 4.7-5.4 covered by our sample, the full band cor- 
responds to the w 2.8-51 keV rest-frame band. All but one of 
the sources observed by Chandra were detected. The Chan- 
dra positions of the detected quasars were found to lie within 
0.2"-0.8" of their optical positions (see Table 1); this is consis- 
tent with the expected positional error. Given the X-ray weak- 
ness of SDSS 0756+4502 (two full-band counts), we will dis- 
cuss the reliability of its detection in more detail in §2.1.1. The 
Chandra non-detection, SDSS 1737+5828, will be discussed 
further in §4. 1 . 

The WAVDETECT photometric measurements are shown in 
Table 2; they are consistent with independent measurements 
obtained using manual aperture photometry and a 2" radius cir- 
cular cell. Background is typically negligible in the source ex- 
traction region, except for SDSS 0941+5947, where the obser- 
vation is partially affected by a background flare; this is taken 
into account in Table 2. The only object with enough counts 
to derive a useful power-law photon index via direct X-ray 
spectral analysis is SDSS 0231—0728. Assuming only Galac- 
tic absorption, it has Y = 2.4+Qg using the unbinned data and 
the C-statistic (Cash 1979) with xspec Version 1 1 .2.0 (Arnaud 
1996); errors are at the 90% confidence level for one interest- 
ing parameter (AC = 2.71; Avni 1976; Cash 1979). 7 This pho- 
ton index is consistent, within the errors, with those of z ~ 0-3 
quasars in the rest-frame 2-10 keV band (r « 1.7-2.3; e.g., 
George et al. 2000; Mineo et al. 2000; Reeves & Turner 2000; 
Page et al. 2003). 

The soft-band images of the detected quasars (full-band im- 
age for BR 0305-4957), adaptively smoothed using the algo- 
rithm of Ebeling, White, & Rangarajan (2003), are shown in 
Figure 2. 

2.1.1. Chandra Detection of SDSS 0756+4502 

Similarly to the z = 5.27 quasar SDSS 1208+0010 presented 
in §2.2 of VOla, SDSS 0756+4502 (z = 4.80) shows only two 
X-ray photons (both in the observed 0.5-2 keV band) within 
0.7" of the optical position of the source. The source is de- 
tected by WAVDETECT using a false-positive threshold of 10~ 4 . 
To investigate further whether the SDSS 0756+4502 X-ray de- 
tection can be considered real, we extracted a 400x400 pixel 2 
(« 200x200 arcsec 2 ) region centered on the optical position 
of the quasar, excluding the immediate vicinity of the quasar 
itself. We covered this region with 30,000 circles of 1" ra- 
dius whose centers were randomly chosen. The choice of the 
1" radius aperture appears reasonable, since the full-band 80% 
encircled-energy radius for on-axis sources is w 0.69" (see 
Fig. 6.3 of the Chandra Proposers' Observatory Guide) and 
the average positional error for the other quasars detected by 
Chandra is w 0.5". The counts obtained in each circle were 
averaged, giving 0.060 counts circle -1 in the full band and 

caused by molecular contamination of the ACIS filters, we have applied a time- 
ry response files which have been used in the X-ray spectral analyses presented 
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TABLE 1 

Chandra AND XMM-Newton OBSERVATION LOG 



Object 




Optical 


Optical 


Aopt-x d 


X 


-ray 




Exp. Time " 


HET 






Name 


z 


«2000 


§2000 


(arcsec) 


Obs 


. Date 


(ks) 


Obs. Date 


Ref. 


SDSSp J020651.37+121624.4 


4.81 


02 06 51.4 


+ 12 16 24.4 


0.8 


2001 


Dec 


29 


5.88 


2001 Dec 


10 


(1) 


SDSSp J023 137.65-072854.5 


5.41 


02 31 37.6 


-07 28 54.5 


0.2 


2002 


Sep 


27 


4.20 






(1) 


BRB0305-4957 


4.73 


03 07 22.9 


-49 45 47.8 


0.4 


2002 


Sep 


22 


3.71 






(2,3) 


SDSSp J033829. 31+002156. 3 c 


5.00 


03 38 29.3 


+00 21 56.3 




2002 


Feb 


22 


5.49 






(4) 


RD 657 c 


4.96 


03 38 30.0 


+00 18 39.9 


1.7 


2002 


Feb 


22 


4.85 






(5) 


SDSSp J075618.14+410408.6 


5.09 


07 56 18.1 


+41 04 08.6 


0.2 


2002 


Feb 


08 


7.33 


2001 Dec 


09 


(1) 


SDSSp J075652.07+450258.9 


4.80 


07 56 52.1 


+45 02 58.9 


0.5 


2002 


Oct 


04 


6.97 


2001 Dec 


10 


(1) 


SDSSp J091316.56+591921.5 


5.11 


09 13 16.6 


+59 19 21.5 


0.7 


2002 


Mar 


07 


9.84 


2002 Feb 


07 


(1) 


SDSSp J094108.36+594725.8 


4.82 


09 41 08.4 


+59 47 25.8 


0.2 


2001 


Oct 


25 


4.16 


2001 Nov 16-18 


(1) 


SDSSp J095151. 17+594556.2 


4.86 


09 51 51.2 


+59 45 56.2 


0.3 


2002 


Apr 


14 


5.13 


2002 Feb 


07 


(1) 


SDSSp J102332.08+633508.1 


4.88 


10 23 32.1 


+63 35 08.1 


0.1 


2002 


Oct 


09 


4.70 






(1) 


SDSSp J173744.87+582829.5 


4.94 


17 37 44.9 


+58 28 29.5 




2002 


Aug 


05 


4.62 


2002 Jun 


08 


(1) 


SDSSp J221644.02+001348.3 


4.99 


22 16 44.0 


+00 13 48.3 


0.4 


2002 


Nov 


01 


7.41 


2002 Jun 


14 


(1) 



NOTE. — The optical information for the quasars presented here can be found in the papers cited in the reference column. For BR B0305— 4957 we 
use the redshift z = 4.73 from Peroux et al. (2001) derived from two emission lines (Si TV-0 IV] and C IV) in their high-quality optical spectrum. 
Storrie-Lombardi et al. (2001) reported z = 4.78 (see their Table 1) and z = 4.82 (see their Fig. 3) derived from a lower-quality optical spectrum 
with less wavelength coverage. 

a Distance between the optical and X-ray positions; a blank entry indicates no X-ray detection. b The Chandra exposure time has been corrected for 
detector dead time. The XMM-Newton exposure time, computed from the exposure maps (i.e., taking into account the decreasing effective exposure 
at increasing off-axis angles from the aimpoint), is for the observation after the background flares have been removed (see §2.2). c XMM-Newton 
observation. 

REFERENCES. — Anderson et al. 2001; (2) Peroux et al. 2001; (3) Storrie-Lombardi et al. 2001; (4) Fan et al. 1999; (5) Djorgovski et al. 2003. 

themselves) with those obtained from the ROSAT deep sur- 
vey of the Lockman Hole (Hasinger et al. 1998) and the 
XMM*wton moderately deep survey HELLAS2XMM (Baldi 
et al. 2002) indicates that in the proximity of the SDSS 
quasars there is no evidence for an X-ray source overdensity. 
We found a surface density of N(> 5)=751+|q} deg~ 2 for the 
X-ray sources surrounding our SDSS targets (where S is the 
soft X-ray flux limit we reached with the present observations: 
«2 x 10~ 15 ergcirr 2 s" 1 ). 

2.1.3. X-ray Extension of the Quasars 

To constrain the presence of gravitational lensing (e.g., 
Barkana & Loeb 2000; Comerford, Haiman, & Schaye 2002; 
Wyithe & Loeb 2002), X-ray jets close to the quasars (e.g., 
Schwartz 2002), and X-ray scattering by dust in the intergalac- 
tic medium (e.g., Telis et al. 2003), we performed an analy- 
sis of X-ray extension for all the quasars observed by Chandra 
following the procedure described in §3.5 of V03. All of the 
quasars' X-ray images are consistent with the sources being 
point-like. The separation between multiply imaged sources 
at z w 4-5 is expected to be « 1-2" (Barkana & Loeb 2000). 



0.022 counts circle in the soft band. The Poisson probabil- 
ities of obtaining two counts or more when 0.060 and 0.022 
counts are expected are w 1.7 x 10~ 3 and w 2.4 x 10~ 4 , cor- 
responding to one-tailed Gaussian probabilities of ~ 3.3o and 
w 3.9o, respectively. 

We have directly verified the validity of the Poisson ap- 
proximation for the SDSS 0756+4502 field using the count 
distribution from the 30,000-circle analysis described above. 
The results of our simulations are in reasonable agreement 
with the Poisson probabilities. Therefore, in the following, 
SDSS 0756+4502 will be considered detected by Chandra. 

2.1.2. X-ray Sources Close to the Quasars 

Following the method used in §3.4 of V03, we have 
searched for possible companions for all of the objects ob- 
served by Chandra over a region of xs 100"xl00" centered 
on the quasars (see Fig. 2), corresponding to a physical scale 
of w 630 x 630 kpc 2 at the average redshift of our sample 
(z « 5.0). Approximately one faint serendipitous X-ray source 
is found in each quasar field. A comparison of their 0.5-2 keV 
cumulative number counts (excluding, of course, the quasars 



TABLE 2 

Chandra SOURCES: X-RAY COUNTS 







X-ray Counts" 




Object 


[0.3-0.5 keV] 


[0.5-2 keV] 


[2-8 keV] 


[0.5-8 keV] 


SDSS 0206+1216 


<3.0 


3.y_ 24 


<4.8 


6 9+ 18 


SDSS 0231-0728 




22.0+58 




23.9+|° 


BR 0305-4957 


<3.0 


2.0+fJ 


2.<£?;? 


4.0+t 2 . 


SDSS 0756+4104 


2 0+ 2 ' 7 
z.u_, 3 


10.0+ 4 / 3 


3.9I 3 ; 2 


14.8+ 4 ^ 


SDSS 0756+4502 


<4.8 


2.0+!? 


< 3.0 


2.0+fl 


SDSS 0913+5919 


<4.8 


s-atii 


<4.8 


6 0+ 36 


SDSS 0941+5947 


<3.0 


4.8+J;? 


<3.0 


4 o+ 3 ' 2 
^■ z -2.0 


SDSS 0951+5945 


<3.0 


3.oi 2 ;? 


<4.8 


t q+3.2 


SDSS 1023+6335 


<3.0 


3 +2 9 


<4.8 


4 0+ 3 ' 2 


SDSS 1737+5828 


<3.0 


<3.0 


<3.0 


<3.0 


SDSS 2216+0013 


<4.8 


6 0+ 3 6 
o.u_ 24 


<4.8 


6 n+3.8 



" Errors on the X-ray counts were computed according to Tables 1 and 2 of Gehrels 
(1986) and correspond to the la level; these were calculated using Poisson statistics. 
The upper limits are at the 95% confidence level and were computed according to 
Kraft, Burrows, & Nousek (1991). 
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FIG. 2.— Chandra 0.5-2 keV (0.5-8 keV for BR 0305-4957) and XMM-Newton 0.3-12 keV (for RD 657) images of the X-ray detected quasars presented in 
this paper. In each panel, the horizontal axis shows the Right Ascension, and the vertical axis shows the Declination (both in J2000 coordinates). North is up, and 
East to the left. The images are 97" x 97" and have been adaptively smoothed at the 2a level [except for BR 0305-4957, RD 657, SDSS 0951+5945 (1.5o level), 
and SDSS 0756+4502 (1.3a level)] using the algorithm of Ebeling et al. (2003). Crosses mark the optical positions of the quasars. The image of RD 657 does not 
cover the position of SDSS 0338+0021, since the two quasars are separated by 196". 



Typical flux ratios of multiply imaged quasars at high redshift strongly dependent on the shape, mass, redshift, and location of 
are w 2-8 (G. Chartas 2003, private communication), being the lensing system. Given the limited number of counts for our 



X-RAY EMISSION FROM z > 4.7 QUASARS 



5 



detected sources, we expect that only in w 30% of the present 
sample would we be able to detect such multiple images. 

2.1.4. X-ray Variability of the Quasars 

Given tentative reports of increasing quasar X-ray variabil- 
ity with redshift (in the sense that quasars of the same X-ray 
luminosity are more variable at z > 2; Manners, Almaini, & 
Lawrence 2002), we have searched for X-ray variability by ana- 
lyzing the photon arrival times in the 0.5-8 keV band using the 
Kolmogorov-Smirnov test. No significant variability was de- 
tected from any of our quasars. We note, however, that the anal- 
ysis of X-ray variability for the quasars studied here is limited 
by both the short X-ray exposures in the rest frame (« 15 min; 
see Table 1) and the low X-ray fluxes of most of our quasars. 

2.2. XMM-Newton Observation 

The SDSS 0338+0021 field was observed by XMM-Newton 
for a nominal exposure time of 20.0 ks with the European Pho- 
ton Imaging Camera (EPIC) pn (Striider et al. 2001) and for 
22.3 ks with the EPIC MOS (Turner et al. 2001) instruments. 
The data obtained with the Reflection Grating Spectrometer 
(RGS; den Herder et al. 2001), with a nominal exposure time of 
22.8 ks, will not be considered in the following given the weak- 
ness of the target. Both the pn and MOS data were acquired 
in the "full-frame mode" using the thin filter (due to the likely 
weakness of the target in the X-ray band). The data were re- 
duced with the XMM-Newton Science Analysis Software (SAS, 
Version 5.3.3) with the latest calibration products; we used the 
EPCHAIN and EMPROC tasks to generate valid photon lists. We 
filtered the data to avoid background flares (due to « 100 keV 
protons; e.g., De Luca & Molendi 2002) which affect most of 
this XMM-Afewton observation. To select good time intervals, 
we apply thresholds of 0.35 counts s _1 and 0.15 counts s _1 to 
pn and MOS data in the 10-12.4 keV and 10-13 keV energy 
bands, respectively (A. Baldi 2002, private communication). 
In these energy bands, the counts are mostly from the back- 
ground. The final, on-axis useful exposure time for this obser- 
vation is 5.5 ks forpn and MOS. Both single-pixel and double- 
pixel events (patterns 0—4) were used when extracting the pn 
counts, while patterns 0-12 were used for the MOS. 

We performed source detection in the 0.3-2 keV, 0.5-2 keV, 
0.5-10 keV, 2-10 keV, and 0.3-12 keV bands using the tasks 
EBOXDETECT and EMLDETECT and adopting a detection like- 
lihood threshold of 10, corresponding to the « 4a detection 
level. We did not consider the energy range below 0.3 keV be- 
cause of the enhanced background at low photon energies. The 
target, SDSS 0338+0021, was detected in none of the bands; 
its X-ray non-detection will be discussed further in §4. 1 . 

After our XMM-Newton observation was assigned, Djorgov- 
ski et al. (2003) found az = 4.96 quasar (RD 657) located 
196" from the optical position of SDSS 0338+0021. In the 
XMM-Newton pn field-of-view we found a 4.3a X-ray source 
detected 1.7" from the optical position of RD 657 (see Table 1) 
in the 0.3-12 keV band (15.3±^ counts). We note the XMM- 
Newton positional error is larger (« 4-5"; Metcalfe 2002) than 
that of Chandra; therefore the X-ray and optical positions are 
consistent. Assuming the integral 0.5-2 keV number counts 
of Hasinger et al. (1998) and an X-ray error circle of 1.7", 
we expect s=a 2.5 x 10~ 4 spurious associations given the soft 

8 HET observations were not performed for SDSS 0231-0728, BR 0305-4957, 

9 The only exception is SDSS 0756+4502, observed by the HET 10 months befc 



X-ray flux of this source (3.6xl0~ 15 ergcm~ 2 s _1 ; see Ta- 
ble 3). Therefore, we are confident that the serendipitously 
detected X-ray emission is associated with RD 657. For this 
source we compared the 4.3a detection obtained by the SAS 
detection tasks with that obtained in a different way with S EX- 
TRACTOR (Bertin & Arnouts 1996), finding similar results in 
both the number of counts and the significance of the source 
(3.6a). 

Inspection of the X-ray image at the source position con- 
firms the X-ray detection (see Fig. 2 for the 1.5a level adap- 
tively smoothed image in the 0.3-12 keV band). RD 657 also 
seems to be present in the 0.3-2 keV adaptively smoothed im- 
age, although the quasar is not detected in this band by the 
SAS procedure described above. Even though there has been 
speculation that SDSS 0338+0021 and RD 657 mark a large- 
scale structure, perhaps a proto-cluster at z ~ 5 (Djorgovski 
et al. 2003), we found no clear evidence for diffuse X-ray 
emission or an over-density of X-ray sources around these two 
quasars. Our constraints are not tight, however, due to the 
short exposure time and relatively high background of XMM- 
Newton. Since the background in the XMM-Newton obser- 
vation is higher and more complex than in the Chandra ob- 
servations, for SDSS 0338+0021 and RD 657 we do not re- 
port the counts in individual bands but directly derived the 
main X-ray parameters in Table 3 (see §3). We have also 
analyzed the XMM-Newton data with less conservative back- 
ground screening choices and our main results are unchanged; 
SDSS 0338+0021 is not detected and RD 657 is detected. 

2.3. Hobby-Eberly Telescope Observations 

Given that quasars frequently vary in both the optical and 
X-ray bands, one must be concerned about the reliability 
of the frequently employed optical-to-X-ray power-law slope 
whenever observations in the two bands are not simultaneous. 
To search for optical variability, eight of the 13 quasars in 
our sample were observed with the 9-m Hobby-Eberly Tele- 
scope (Ramsey et al. 1998). 8 Short (2 min; 4 min for 
SDSS 0941+5947) /-band images of the quasar fields were 
taken with the HET's Marcario Low-Resolution Spectrograph 
(LRS; Hill et al. 1998a,b; Cobos Duenas et al. 1998) generally 
within 1-2 months of the X-ray observations (1-2 weeks in the 
rest frame; see Table l). 9 The HET images covered fields « 4 
on a side, and using the published finding charts of the fields it 
was possible to obtain an approximate photometric calibration 
using stellar objects in the fields. Although it is difficult to give 
precise estimates given the difference between the LRS and 
SDSS bandpasses, it is clear that none of the quasars observed 
with the HET has varied significantly (by ^ 30%) in brightness 
from its published value. These constraints on the optical flux 
variability allow for at most an uncertainty of « 0.05 in the a ox 
determination (see §3). The lack of large brightness variations 
in the sample is not surprising considering the relatively small 
amount of time in the rest frame (several months) between the 
SDSS photometry and the X-ray/HET observations. 

To investigate the nature of SDSS 1737+5828, the only 
quasar that was not detected in the Chandra observations, we 
performed a 30-min LRS spectroscopic observation. This was 
taken with the 300 line mm -1 grating, an OG515 blocking fil- 
ter, and a slit width of 2" (with a seeing of « 1 .6"), providing a 
spectrum in the wavelength range 5100-10200 A at a resolving 

DSS 0338+0021, RD 657, and SDSS 1023+6335. 
the X-ray observation. 
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power of w 240-300. The spectroscopic data were reduced us- 
ing standard IRAF routines. The wavelength calibration was 
performed using a Ne arc lamp, while the flux calibration was 
achieved using the standard star BD +17 4708, observed dur- 
ing the same night. The HET spectrum (black curve in Fig. 3) is 
discussed in §4.1 along with the SDSS spectrum (adapted from 
Anderson et al. 2001; grey curve in Fig. 3). 



, , , | , , , , | r 

Lya a fas 




Wavelength (A) 



FIG. 3.— Optical spectra of SDSS 1737+5828 obtained with the HET (30- 
min exposure; black curve) and by the SDSS (adapted from Anderson et al. 
2001; grey curve). No correction for the telluric lines has been applied to the 
HET spectrum (the affected regions of the HET spectrum are shaded). The 
likely absorption features from N V and Lya are shown (see §4. 1 for details). 

A direct comparison of the continuum levels of the two spec- 
tra is difficult. While the uncertainty in the calibration of the 
SDSS spectrum is less than w 20%, it is hard to obtain an ac- 
curate HET flux calibration mostly because of the variations of 
effective aperture produced by its optical design (in addition to 
problems due to slit losses). 

3. MULTI-WAVELENGTH PROPERTIES OF Z > 4.7 QUASARS 

The principal optical, X-ray, and radio properties of the 
quasars in our sample are shown in Table 3. A description is 
as follows: 

Column (1 ). — The abbreviated name of the source. 

Column (2). — The Galactic column density (from Dickey & 

Lockman 1990) in units of 10 20 cm~ 2 . 

Column (3). — The monochromatic rest-frame AB 1450 (i +z ) 
magnitude (defined in §3b of Schneider et al. 1989) with es- 
timated errors w 0. 1 mag. The absolute photometry is correct 
to within w 0.05 mag (Anderson et al. 2001), while the compu- 
tation of the AZ?!45o(i +z ) magnitude using the photometric mea- 
surements and the composite SDSS quasar spectrum (Vanden 
Berk et al. 2001) is correct to within « 0.05-0.10 mag. 
Columns (4) and (5). — The 2500 A rest-frame flux density 
and luminosity. These were computed from the AB U50 n +z \ 
magnitude assuming an optical power-law slope of a = —0.79 
(5 V °= V 01 ; Fan et al. 2001) to allow direct comparison with the 
results presented in VOla and V03. The 2500 A rest-frame 
flux densities and luminosities are reduced by k 15% when 



the power-law slope of the optical continuum is changed to 
a = -0.5 (e.g., Vanden Berk etal. 2001; Schneider et al. 2001). 
Column (6). — The absolute B-band magnitude computed us- 
ing a = —0.79. If a = —0.5 is adopted for the extrapolation, 
the absolute B-band magnitudes are fainter by « 0.35 mag. 
Columns (7) and (8). — The observed count rate in the 
0.5-2 keV band and the corresponding flux (Fx), corrected 
for Galactic absorption. This flux has been calculated using 
PIMMS (Version 3. 2d; Mukai 2002) and a power-law model 
with r = 2.0, as derived for samples of z ~ 0-3 quasars (e.g., 
George et al. 2000; Mineo et al. 2000; Reeves & Turner 2000; 
Page et al. 2003) and z > 4 optically luminous, PSS quasars 
(V03); see also §4.3. 

The X-ray fluxes reported in this paper have been corrected 
for the ACIS quantum efficiency decay at low energy (see 
Footnote 7). For each quasar, we produced an ancillary re- 
sponse file (ARF) assuming an exponential decay for the ACIS 
quantum efficiency. From each ARF we derived an effec- 
tive area that was used to replace the original ACIS-S effec- 
tive area in PIMMS. The 0.5-2 keV flux correction is w 20% 
for the Chandra observations presented here. All of the pa- 
rameters derived from the measured soft X-ray flux have been 
corrected. Following this procedure, we also corrected all 
of the X-ray fluxes of the z > 4 quasars previously observed 
by Chandra (see Table Al in the Appendix for an updated 
version of the relevant numbers for z > 4 quasars; also see 
http://www.astro.psu.edu/users/niel/papers/). 

The soft X-ray flux changes by w 1-2% over a reasonable 
range in the adopted power-law photon index, T = 1 .7-2.3. The 
soft X-ray flux derived from the 0.5-2 keV counts is generally 
similar (to within w 5-20%) to that derived from the 0.5-8 keV 
counts. 

For RD 657, the 0.5-2 keV flux and subsequent X-ray pa- 
rameters reported below have been derived using the counts in 
the 0.3-12 keV band, where it was detected. 
Columns (9) and (10). — The rest-frame 2 keV flux density 
and luminosity, computed assuming F = 2.0. 
Column (11). — The 2-10 keV rest-frame luminosity, cor- 
rected for Galactic absorption. 

Column (12). — The optical-to-X-ray power-law slope, a ox , 
defined as 

„ _ log(/2kev// 2 500 A) m 
log(V 2 kev/V 2 500 A) 

where fa k e v and / 2 500 A are tne rest-frame flux densities at 
2 keV and 2500 A, respectively. The « la errors on a ox have 
been computed following the "numerical method" described in 
§ 1.7.3 of Lyons (1991). Both the statistical uncertainties on 
the X-ray count rates and the effects of possible changes in the 
X-ray (r w 1.7-2.3) and optical (a w -0.5 to -0.9; Schneider 
et al. 2001) continuum shapes have been taken into account (see 
§3 of VOla for further details). Changing the power-law slope 
of the optical continuum from a = —0.79 to a = —0.5 induces 
a small increase in the a ox values (more positive by w 0.03). 
Column (13). — The radio-loudness parameter (e.g., Keller- 
mann et al. 1989), defined as R=fs GHz/ /4400 A ( rest frame). 
The rest-frame 5 GHz flux density was computed from the 
FIRST (Becker, White, & Helfand 1995) or NVSS (Condon 
et al. 1998) observed 1.4 GHz flux density assuming a radio 
power-law slope of a = —0.8. The upper limits reported in 
the table are at the 3a level. The rest-frame 4400 A flux den- 
sity was computed from the AB 14 5 ( 1+Z ) magnitude assuming 
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Note. — Luminosities are computed using Hq=70 km s _1 Mpc~', Hm=0.3, and I2 A =0.7. 
"From Dickey & Lockman (1990) in units of 10 20 cm~ 2 . 
b Rest-frame 2500 A flux density in units of 10~ 27 erg cm~ 2 s _1 Hz -1 . 
c Observed count rate computed in the 0.5-2 keV band in units of 10~ 3 counts s _1 . 

d Galactic absorption-corrected flux in the observed 0.5-2 keV band in units of 10~ 15 erg cm~ 2 s -l. These fluxes and the following X-ray parameters have been corrected for the 
ACIS quantum efficiency decay at low energy. 

e Rest-frame 2 keV flux density in units of 10~ 32 erg cm~ 2 s _1 Hz -1 . 

f Rest-frame 2-10 keV luminosity corrected for Galactic absorption in units of erg s _1 . 

g Errors have been computed following the "numerical method" described in § 1.7.3 of Lyons (1991); both the statistical uncertainties on the X-ray count rates and the effects of the 
observed ranges of the X-ray and optical continuum shapes have been taken into account (see the text for details; see also §3 in VOla). 

h Radio-loudness parameter, defined as R = f s GHz// 4400 a ( rest f rame ; e -g-> Kellermann et al. 1989). The rest-frame 5 GHz flux density is computed from the observed 1.4 GHz flux 
density assuming a radio power-law slope of a = -0.8 with / v « v a . 

; 1.4 GHz flux density from NVSS (Condon et al. 1998). 

> 1 .4 GHz flux density from FIRST (Becker et al. 1995). 

k No radio data are available. 

'The soft X-ray flux and the subsequent X-ray parameters have been computed from the 0.3-12 keV counts, where the source is detected. 
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an optical power-law slope of a = —0.79. Radio-loud quasars 
(RLQs) have R > 100, whereas radio-quiet quasars (RQQs) are 
characterized by R < 10 (e.g., Kellermann et al. 1989). 

SDSS 0913+5919 is the only radio-loud object (R w 41 1) in 
our sample and is also the highest redshift (z = 5.11) RLQ thus 
far detected in the X-ray band. For one quasar, BR 0305—4957, 
no radio data are available. None of the remaining quasars is a 
RLQ, although in some cases the observed-frame 1 .4 GHz flux- 
density constraints provided by the NVSS are not sufficiently 
tight to prove that a quasar is radio quiet (R < 10); at most 
these objects could be radio moderate (see Table 3). Given that 
« 85-90% of the quasar population is radio quiet (e.g., Stern 
et al. 2000), we expect the majority of the quasars with radio- 
loudness upper limits of 10-100 to be radio quiet. 

4. DISCUSSION 

Thanks mostly to Chandra (e.g., VOla; Brandt et al. 2002; 
Bechtold et al. 2003; V03) and ROSAT (e.g., KBS00; V03) 
observations, we have seen a substantial growth in the data 
needed to study the X-ray properties of z > 4 quasars. At 
present, there are w 70 quasars at z > 4 with X-ray detections. 10 
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FIG. 4. — Observed-frame, Galactic absorption-corrected 0.5-2 keV flux 
versus ABi4 5 o(i +z ) magnitude for z > 4 AGNs. Object types are shown in the 
key. In particular, the objects presented in this paper are plotted as filled trian- 
gles {Chandra detections) and a filled square (XMM-Newton detection), while 
large, thick, downward-pointing arrows indicate X-ray upper limits for the ob- 
jects of our sample. Chandra RLQs are plotted as circled triangles. ROSAT 
upper limits (small downward-pointing arrows) are at the 3o confidence level. 
The slanted dashed lines show the a ox =— 1 .5 and a ox =— 1.8 loci at z = 5.0 (the 
average redshift of the present sample). The dotted line shows the best-fit cor- 
relation reported in §4 [Equation (2)] for the whole sample of z > 4 optically 
selected RQQs. 

Figure 4 shows the observed-frame, Galactic absorption- 
corrected 0.5-2 keV flux versus AB 145 o(i_|_ z ) magnitude for a 
compilation of z > 4 Active Galactic Nuclei (AGNs). The de- 
tected objects presented in this paper are shown as filled trian- 
gles (10 Chandra detections) and a filled square (one XMM- 
Newton detection), while large and thick downward-pointing 
arrows indicate the two quasars in our sample with X-ray up- 
per limits. Because this figure is constructed from directly ob- 



served physical quantities, it is robust and useful for planning 
future X-ray observations of high-redshift quasars. We present 
a complementary discussion based on a ox , a parameter some- 
what less directly tied to observation, in §4.2. 

Quasars are well-known X-ray emitters at z ~ 0-3 (e.g., Vi- 
gnali et al. 1999; Reeves & Turner 2000; Page et al. 2003). 
Figure 4 indicates that X-ray emission is also a universal prop- 
erty of z > 4 quasars. 

Previous studies of z > 4 quasars have shown that there is 
a significant correlation between AB 1450 ( 1+Z ) and 0.5-2 keV 
flux (VOla; V03). The new X-ray detections and upper limits 
provided by the present work significantly improve the study 
of this correlation at fainter optical magnitudes. The average 
ABi45o(i +z ) magnitude of the sample presented here is 19.7, 
whereas the average for the bright PSS sample of V03 is 17.9. 
The presence of this correlation indicates that the same engine 
(namely accretion onto a supermassive black hole) is powering 
both the ultraviolet and X-ray emission. To study this corre- 
lation using all the available X-ray detections (41) and upper 
limits (52) for z > 4 optically selected RQQs (the majority of 
the objects plotted in Fig. 4), we have used the ASURV soft- 
ware package Rev 1.2 (La Valley, Isobe, & Feigelson 1992). 
To evaluate the significance of the correlation, we used sev- 
eral methods available in ASURV, namely the Cox regression 
proportional hazard (Cox 1972), the generalized Kendall's T 
(Brown, Hollander, & Korwar 1974), and the Spearman's p 
models. We have chosen to use only optically selected RQQs 
in order to have a homogeneous sample (see §4 of Vignali et 
al. 2002, hereafter V02, for the results of a sample of X-ray 
selected quasars; also see Silverman et al. 2002 and Castander 
et al. 2003). Once the Broad Absorption Line quasars (BALQ- 
SOs), which are known to be absorbed in the X-ray band (e.g., 
Brandt, Laor, & Wills 2000; Gallagher et al. 2001, 2002; Green 
et al. 2001; VOla), are excluded from the analysis (there are 
four in the original sample), we find that the AB 1450 (i +z ) versus 
Fx correlation is significant at the > 99.9% level (all the tests 
above provide similar results). According to the EM (Estimate 
and Maximize) regression algorithm (Dempster, Laird, & Ru- 
bin 1977), the correlation is parameterized by 

log(F x /ergcm- 2 s- 1 ) = [-(0.33±0.05)AB 1450(1+ , ) -(7.89 ±0.91)], 

(2) 

which is plotted as a dotted line in Figure 4. Since 
AB 14 5o(i +z )=— 2.5 log / 1450 ^—48.6, the previous equation may 
also be written as 



logOFx/ergcm^s- 1 ) - (0.83±0.13) log/ 1450 



(3) 



In Figure 5 the quantity /*x,obs/^x,fit is shown as a func- 
tion of AB 1450 (i +z ) magnitude for the z > 4 optically se- 
lected RQQs observed by Chandra. Fx.obs represents the 
measured, Galactic absorption-corrected 0.5-2 keV flux, and 
Fxfit the flux expected from the best-fit correlation above. 
The majority of the z > 4 quasars observed by Chandra are 
within a factor of two of the best-fit correlation. The only 
quasar clearly showing a large difference between the mea- 
sured and expected fluxes is SDSS 0231—0728 (the bright- 
est X-ray source in our sample). Since its X-ray emis- 
sion is unlikely to be contaminated by nearby X-ray sources 
(see Fig. 2), it appears plausible that its unusually large 
X-ray flux is due to an episode of variability. Unfortu- 
nately, the lack of an HET observation for SDSS 0231-0728 



1 See http://www.astro.psu.edu/users/niel/papers/highz-xray-detected.dat for a regularly updated compilation of X-ray detections at z > 4. 
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does not allow us to examine this hypothesis further (as- 
suming that the X-ray and optical variability are correlated). 
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FIG. 5. — Plot of Fx.obJFxM as a function of AB i450 ( l+z) magnitude for 
the z > 4 optically selected RQQs observed by Chandra. Fx.obs represents the 
Galactic absorption-corrected 0.5-2 keV flux observed by Chandra and F x ,flt 
the flux expected on the basis of the best-fit correlation between the soft X-ray 
flux and AB 1450 ( 1+Z ) magnitude discussed in §4 [Equation (2)]. Symbols are 
the same as in Figure 4. The dotted line indicates where the measured and 
expected fluxes are the same. To avoid symbol crowding, we have only shown 
a representative average uncertainty on the flux ratios (the error bar near the 
middle right). 



A constant fit to the data points of Figure 5 (including the 
associated errors, but excluding the upper limits) is statistically 
unacceptable (% 2 = 54.6 for 33 degrees of freedom). Excluding 
SDSS 0231-0728 improves the quality of the fit (% 2 = 45.3), 
although it remains unacceptable. This indicates that there is 
a real dispersion in the quantity plotted in Figure 5 above the 
statistical noise (see also §4.2). 

4.1. X-ray Non-Detections: Likely Explanations 

All but two (SDSS 0338+0021 and SDSS 1737+5828) of 
the quasars presented here have been detected in the X-ray 
band. The non-detection of SDSS 0338+0021 is probably 
due to the short exposure time of the XMM-Newton observa- 
tion once the flaring-background intervals have been removed. 
Assuming the best-fit parameters of the ABn^i +7 ^ versus 
0.5-2 keV flux correlation reported in Equation (2), we would 
expect a 0.5-2 keV flux of w 3.2 x 10~ 15 erg cm" 2 s _1 . Al- 
though SDSS 0338+0021 was observed at the aimpoint, the on- 
axis flux limit of the XMM-Newton EPIC pn in a w 5.5 ks ex- 
posure (see Table 1) is a factor of « 2 higher than that achieved 
with Chandra in the same amount of time. 

No Chandra source counts were detected at the optical po- 
sition of SDSS 1737+5828. The faint ROSAT source reported 
as a possible counterpart of SDSS 1737+5828 by Anderson et 
al. (2001) is clearly present in the Chandra observation, but its 
X-ray position is 45" from the quasar's optical position. There 
are two possible explanations for the SDSS 1737+5828 non- 
detection: this quasar is either strongly absorbed or intrinsi- 

11 This average value is flatter by 0.04 than that originally reported in V03 because 



cally X-ray weak. If the non-detection is due to the presence 
of intrinsic X-ray absorption, a column density of Nh ~ 1.2 x 
10 23 cm~ 2 is necessary to reproduce the Chandra constraints 
(see Tables 2 and 3) assuming the average a ox =— 1 -71 of the 
sample (see §4.2) and an X-ray photon index of F = 2.0. The 
absorption possibility finds partial support from the published 
SDSS 1737+5828 optical spectrum (Anderson et al. 2001; see 
the grey curve in Fig. 3), which shows evidence of a strong, 
narrow N V absorption feature at z « 4.92. This absorption 
feature has a fairly high formal significance in the SDSS spec- 
trum, but given the importance of the intrinsic absorption inter- 
pretation, we wished to acquire an independent confirmation of 
the line; the new observation would also allow an investigation 
of any line variability. We therefore obtained a 30-min spec- 
troscopic observation of SDSS 1737+5828 with the HET (see 
§2.3). The resulting spectrum (black curve in Fig. 3) confirms 
the presence of a clear narrow N V absorption feature. This 
line is unresolved (the spectral resolution is w 1000 km s -1 ) 
and has a blueshift of w 700 km s -1 relative to the quasar red- 
shift of 4.94 reported by Anderson et al. (2001). In both the 
SDSS and HET spectra we see evidence for narrow Lya absorp- 
tion spanning the redshift of the N V absorption feature with 
the strongest absorption appearing at z ~ 4.90 and z ~ 4.945. 
There may be features at similar redshifts in Si IV and C IV, 
but the signal-to-noise ratio of the spectra prevents us from un- 
ambiguously identifying these systems. Nevertheless, the N V 
line alone strongly suggests that significant internal absorption 
is present in SDSS 1737+5828, since N V absorption of this 
strength is rarely caused by intervening galaxies (e.g., York et 
al. 1991) and the absorption is blueshifted with respect to the 
emission-line redshift. The Lya emission line appears weak; 
the rest-frame equivalent width is only w 35 A. On average, 
the rest-frame equivalent width of Lya + N V in high-redshift 
quasars is « 70-80 A (Schneider, Schmidt, & Gunn 1991; Fan 
et al. 2001). The evidence from the HET and SDSS spectra 
strongly suggests that significant internal absorption is present 
in SDSS 1737+5828, but the quasar is clearly not a strong high- 
redshift BALQSO such as SDSS 0856+5252 (z = 4.79; Ander- 
son et al. 2001), SDSS 1044-0125 (z = 5.74; Fan et al. 2000; 
Brandt et al. 2001; Maiolino et al. 2001; Goodrich et al. 2001), 
or SDSS 1 129-0142 (z = 4.85; Zheng et al. 2000). 



4.2. a ox Results 

Previous reports (e.g., VOla; Bechtold et al. 2003; V03) 
of steeper a ox indices for z £ 4 quasars than for lower- 
luminosity, lower-redshift quasars are confirmed by the present 
analysis; see Figure 6. The 12 RQQs of our sample have 
(oc ox )=— 1.71+0.05 (the error represents the standard deviation 
of the mean). This average a ox value is considerably steeper 
than that obtained for, e.g., the Bright Quasar Survey (BQS; 
Schmidt & Green 1983) RQQs at z < 0.5 ((a ox }=- 1.56+0.02; 
see Brandt et al. 2000 and VOla). The sample of eight optically 
luminous z > 4 PSS RQQs presented in V03 has an even steeper 
average a ox ((a ox )=-1.77±0.03). u The PSS sample of V03 
has (log (/2500 k/ er S s_1 Hz -1 )) = 32.2, while the one pre- 
sented here has (log (/ 2 soo A/ er 8 s_1 Hz -1 )} = 31.5. Accord- 
ing to the known anti-correlation between a ox and 2500 A lu- 
minosity (e.g., Avni, Worrall, & Morgan 1995; Vignali, Brandt, 
& Schneider 2003, hereafter VBS03; also see the discussion 

of the flux correction discussed in §3. 
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later in this subsection) following the equation 12 

a ox = -(0.095 ±0.021) x / 250oA + (1.295 ±0.648), (4) 

we expect an average a ox of —1.70 for the sample of RQQs 
presented here, in good agreement with the average a ox value 
reported above. A summary of these results is shown in Fig- 
ure 6: the a ox distribution for the BQS RQQs (top panel) is 
clearly different from that of the z > 4 RQQs observed by 
Chandra thus far (bottom panel). In both a ox distributions, 
some of the quasars with the steepest a ox values are BALQSOs 
or mini-BALQSOs; areas with arrows indicate upper limits. 



2003; V03; Priddey et al., in preparation); only one of these is 
not detected in the X-ray band (VOla). 

TABLE 4 

Principal Properties of Three Optically Selected Samples 
of rqqs at low and high redsh1fts 





Of 



23 



BQS RQQs 



Chandra RQQs al z>4 



v///x$m%\ , : 



Sample Name 


N 

^detections 
^upper limits 


BQS 
46 
45 
1 


SDSS 
16 
13 

3 


PSS±BRI 
19 
18 
1 


Range of z 
Median z 


0.061-0.472 
0.17±0.01 
0.16 


4.62-6.28 
5.11±0.11 
4.96 


4.00^173 
4.28±0.04 
4.23 


Range of log / 

, 6 6 2500 A 

log / 
* & 2500 A ' 

Median log / 2500 /{ 


29.14-30.70 


31.22-31.75 


31.54-32.59 


29.66±0.05 


31.53±0.04 


31.96±0.06 


29.60 


31.55 


31.92 


Range of log L 2 -iokeV 

< log 7-2-10 keV > " 

Median log L 2 _ 10 kcV 


43.36-44.79 
44.07±0.10 
44.13 


44.53^15.89 
44.99±0.09 
45.03 


44.75^15.82 
45.34±0.07 
45.41 


Range of oc ox 
{ a,» > " 
Median a ux 


-1.23 to -2.06" 
-1.56 ±0.02 
-1.55 


-1.43 to -1.82 
-1.69 ±0.03 
-1.71 


-1.56 to -1.94 
-1.73 ±0.02 
-1.73 


References 


(1,2,3,4) 


(4,5,6,7,8) 


(5,7,9) 



-2.5 



-1.5 



a 



NOTE. — The luminosities are in the quasars' rest frames and have been 
obtained (or converted) adopting Hq=70 km s -1 Mpc~' in a A-cosmology 
with I2m=0.3 and Q. A =0J. ^ is in units of erg s -1 Hz~' and L 2 _io keV 

is in units of erg s~ 1 . Corrections for the Chandra ACIS quantum efficiency 
decay at low energies have been applied in the computation of the X-ray 
luminosities and a ox values (see §3 for details). 

BALQSOs and mini-BALQSOs have been excluded from all of the sam- 
ples used here. For the SDSS and PSS±BRI samples at z > 4, the values 
presented here have been obtained using Chandra and XMM-Newton re- 
sults. For the BQS sample, the rest-frame 2-10 keV luminosities have been 
obtained from a representative subsample of 21 BQS RQQs observed by 
ASCA (George et al. 2000) and BeppoSAX (Mineo et al. 2000) in the ob- 
served Rs 0.5-10 keV and ks 0.1-10 keV bands, respectively. This subsam- 
ple covers approximately the same range in redshift as the original sample 
of 46 BQS RQQs with M B < -23 used in §3.1.1 of Brandt et al. (2001). 
a The errors on average values represent the standard deviation of the mean. 
b ct ox =-2.06 represents the slope of PG 2214±139. See Appendix A of 
Laor & Brandt (2002) for a discussion about possible absorption in this 
quasar. 

REFERENCES. — (1) Brandt et al. 2000; (2) George et al. 2000; (3) Mineo 
et al. 2000; (4) Brandt et al. 2001; (5) VOla; (6) Brandt et al. 2002; (7) 
Bechtold et al. 2003; (8) this paper; (9) V03. 



FIG. 6. — a ox distributions for the BQS RQQs observed by ROSAT (top) 
and z > 4 RQQs observed by Chandra (bottom). The shaded areas indicate 
the BALQSOs or mini-BALQSOs; the areas with arrows indicate 95% confi- 
dence upper limits. The average uncertainty in the a ox measurements from the 
Chandra observations is shown in the bottom panel as the error bar near the 
upper right. 

The principal properties (redshifts, rest-frame 2500 A and 
2-10 keV luminosities, and a ox values) of representative sam- 
ples of RQQs at low and high redshifts are reported in Table 4. 

We find no significant correlation of oc ox with either rest- 
frame 2500 A luminosity or redshift using only the available 
Z > 4 RQQs. However, we can also use the new X-ray detec- 
tions and upper limits presented here to extend the VBS03 study 
of a ox dependence upon 2500 A luminosity and redshift. 13 As 
in VBS03, we mostly use the SDSS RQQs taken from the Early 
Data Release (EDR) quasar catalog (Schneider et al. 2002). In 
particular, at z < 4 we use the 128 EDR RQQs of VBS03 (see 
their §3.1.1) covered by pointed ROSAT PSPC and HRI obser- 
vations, with an overall X-ray detection fraction of « 48%. At 
z > 4 we use all published Chandra (VOla; Brandt et al. 2002; 
Bechtold et al. 2003) and XMM-Newton (Brandt et al. 2001) 
observations of SDSS RQQs in addition to the ones presented 
in this paper. In total we have 20 SDSS RQQ observations at 
z > 4 (including RD 657, observed in one of the fields surveyed 
by the SDSS; see §2.2 and Djorgovski et al. 2003), 16 of which 
have X-ray detections. To improve the statistics at high lumi- 
nosities and redshifts without introducing any evident bias (see 
§3.2 of VBS03), we also include 15 PSS RQQs and seven BRI 
RQQs at z > 4 observed by Chandra (VOla; Bechtold et al. 

12 This is an updated version of Equation (2) of VBS03, taking into account the correction for the quantum efficiency decay of ACIS 

13 See §2.2 of VBS03 for discussion about the method used to derive 2500 A fluxes from the SDSS filters. 



Therefore, the final sample used in the present analysis is com- 
prised of 170 RQQs (« 58% with X-ray detections) over the 
redshift range z~ 0.2-6.3 and the rest-frame 2500 A lumi- 
nosity range log (/2500 k/ er S s_1 Hz -1 ) ss 29.0-32.6. Eleven 
of these quasars are known BALQSOs; results obtained both 
including and excluding BALQSOs are presented below. As 
luminosity is correlated with redshift in flux-limited samples, 
we use a partial correlation analysis technique developed for 
use with censored data (Akritas & Siebert 1996). Partial 
correlation analyses confirm previous findings (e.g., Avni et 
al. 1995; VBS03): a ox depends primarily upon rest-frame 
2500 A luminosity. The significance of this anti-correlation 
is 3.5-4.0O" (slightly higher than in VBS03), depending upon 
whether PSS/BRI quasars and BALQSOs are included in the 
analyses (see Table 5). 

No significant correlation of a ox with redshift is found. 
These results suggest that the mechanism driving quasar broad- 
band emission is similar in the local and early Universe, with 
no evidence for unusual phenomena in high-redshift quasar nu- 
clei such as "trapping radius" effects (e.g., Begelman 1979) or 
accretion-disk instabilities (e.g., Lightman & Eardley 1974). 
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TABLE 5 

Dependencies of Oo X upon redshift and 2500 A luminosity density 

FROM PARTIAL CORRELATION ANALYSIS 



Sample Name 


Number 


Oiox-Z 


ox 6 2500 




of RQQs 




SDSS only, with BALQSOs 


148 


1.2a 


3.7a 


SDSS only, without BALQSOs 


138 


0.9c 


3.5a 


SDSS + PSS, with BALQSOs 


163 


1.3a 


3.9a 


SDSS + PSS, without BALQSOs 


152 


1.2a 


3.9a 


SDSS + PSS + BRI, with BALQSOs 


170 


1.4a 


4.0a 


SDSS + PSS + BRI, without BALQSOs 


159 


1.1a 


4.0a 



Although we did not find a significant correlation of a ox with 
either rest-frame 2500 A luminosity or redshift using only the 
z > 4 RQQs, the spread of a ox for z > 4 RQQs is unlikely to be 
explained solely by measurement errors. To demonstrate this, 
we have applied the likelihood method described in Section IV 
of Maccacaro et al. (1988) to calculate the intrinsic dispersion 
of a ox (accounting for measurement errors). Ignoring upper 
limits, we find an intrinsic dispersion of d = O.O62+Q Q33 (68% 
confidence errors for two interesting parameters; see Maccac- 
aro et al. 1988). Accounting for upper limits would only in- 
crease the intrinsic dispersion (see the bottom panel of Fig. 6). 
For comparison, the "straight" dispersion (not accounting for 
measurement errors) is o s = 0.1 10. 

SDSS 023 1 -0728 and RD 657 have the flattest a ox values in 
our sample. It is possible (see §4) that they experienced vari- 
ability over the period between the optical and X-ray observa- 
tions. An alternative explanation could be the presence of some 
dust obscuring the optical/ultraviolet radiation but not the X-ray 
emission. The optical spectra of these quasars, however, do not 
show any clear evidence supporting this scenario. Radio emis- 
sion can affect the quasar broad-band spectral energy distribu- 
tions producing flatter a ox values. Although significant radio 
emission can be ruled out for SDSS 0231-0728 (see Table 3), 
the radio-loudness upper limit for RD 657 is not sufficiently 
tight to exclude some radio emission from this quasar. 

SDSS 0913+5919, the only radio-loud object in our sam- 
ple and the highest redshift RLQ with an X-ray detection, 
has the highest R parameter (R 411) and the steepest a ox 
(— 1 .67 ±0.11) among the RLQs observed thus far by Chandra 
at z > 4 ( V01 a; V03). X-ray absorption might be an explanation 
for the steep a ox value, although the detection of five of the six 
source counts in the soft band (see Table 2) suggests a different 
interpretation. The anti-correlation between a ox and rest-frame 
2500 A luminosity (e.g., Avni et al. 1995; VBS03) cannot ex- 
plain its steep (for a RLQ) a ox value, since SDSS 0913+5919 
is also the optically faintest z > 4 RLQ observed by Chandra to 
date. 

4.3. Joint X-ray Spectral Fitting of SDSS Quasars at z> 4.8 

To define the average X-ray spectral properties of high- 
redshift SDSS quasars and to make a comparison with those 
obtained for a sample of nine optically luminous PSS quasars 
at z~ 4.1-4.5 (V03), we have selected a sample of 13 SDSS 
quasars at z > 4.8 detected by Chandra (( z ) = 5.23). Their 
total exposure time is 81.6 ks. The quasars in this sample are 
selected to have more than two counts in the full band; this 
choice excludes SDSS 0756+4502 (presented in this paper) 
and SDSS 1208+0010 (VOla). Five of the quasars used in 
this sample have been previously published: the three highest 

If we increase the number of soft-band counts by 10% to take into account crudely the quantum efficiency decay of ACIS at soft energies, we obtain r=1.93±g;^. 
The small differences in the photon indices obtained via joint spectral fitting and band-ratio analysis (with a 10% increment in the number of soft-band counts) are 
probably due to the fact that only the former retains all the spectral information. 



redshift SDSS quasars observed thus far (SDSS 1030+0524, 
SDSS 1306+0356, and SDSS 0836+0054; Brandt et al. 
2002); SDSS 1204-0021 (Bechtold et al. 2003); and 
SDSS 0211-0009 (VOla). Eleven of the 13 quasars are ra- 
dio quiet; the exceptions are SDSS 0836+0054 (radio moderate 
with R w 10) and SDSS 0913+5919 (radio loud with R w 411; 
see Table 3). In the following we present the joint X-ray spec- 
tral analysis both with and without RLQs, since RLQs often 
appear to be characterized by flatter X-ray photon indices (e.g., 
Cappi et al. 1997). The sample does not appear to be bi- 
ased by the presence of a few high signal-to-noise ratio ob- 
jects. For example, the removal of the radio-moderate quasar 
SDSS 0836+0054, the second brightest source in the present 
sample, does not provide significantly different results (see be- 
low). 

Source counts were extracted from 2" radius circular aper- 
tures centered on the X-ray position of each quasar. The back- 
ground was taken from annuli centered on the quasars, avoid- 
ing the presence of nearby X-ray sources. In total we have 120 
counts (in the rest-frame w 2.5-36 keV band) from the source 
extraction regions and five counts from the background regions 
normalized to the size of those used for the sources; half of the 
background counts are from the SDSS 0941+5947 field (see 
§2.1). All of the ARFs used in the spectral analysis have been 
corrected to account for the low-energy quantum efficiency de- 
cay of ACIS (see §2.1). As in §2.1, spectral analysis was car- 
ried out with XSPEC using the background-subtracted data and 
C-statistic (Cash 1979). Since C-statistic fitting is performed 
using the unbinned data, it retains all the spectral information 
and allows one to associate with each quasar its own Galactic 
absorption column density and redshift (for the fitting with in- 
trinsic absorption). 

Joint spectral fitting with a power-law model (leaving the 
normalizations free to vary) and Galactic absorption provides 
a good fit to the data as judged by the small data-to-model 
residuals and checked with 10,000 Monte-Carlo simulations. 
The resulting photon index is T = 1.84^1$. Excluding the 
two non-RQQs from the joint spectral fitting gives a similar 
result (r = I.86+Q37). These photon indices are consistent 
with those obtained for z ~ 0-3 samples of RQQs in the rest- 
frame 2-10 keV band (e.g., George et al. 2000; Mineo et 
al. 2000; Reeves & Turner 2000; Page et al. 2003) and the 
z ~ 4.1^4.5 optically luminous PSS quasars presented in V03 
(r = 1.98 ±0.16). As a consistency check, we compared the 
photon index obtained from joint spectral fitting with that de- 
rived from band-ratio analysis (i.e., the ratio of the 2-8 keV to 
0.5-2 keV counts; we have assumed for the Galactic column 
density the average value for the 13 quasars weighted by their 
number of counts). In this case we obtain T = 1 .85 



+0.22 14 
-0.20- 



We also constrained the presence of neutral intrinsic ab- 
sorption; the joint-fitting technique described above provides 
a counts-weighted average column density. Solar abundances 
have been assumed, although previous optical studies indicate 
that high-redshift quasar nuclei are often characterized by su- 
persolar abundances of heavy elements (e.g., Hamann & Fer- 
land 1999; Constantin et al. 2002; Dietrich et al. 2003). Dou- 
bling the abundances in the fit gives a reduction in the column 
density of a factor of w 2. No evidence for absorption in ex- 
cess of the Galactic column density has been found, as shown 
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by the 68, 90, and 99% confidence regions plotted in Figure 7. 
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FIG. 7. — Confidence regions for the photon index and intrinsic column den- 
sity derived from joint spectral fitting of the 13 SDSS quasars (11 RQQs, two 
non-RQQs) at z > 4.8 thus far detected by Chandra with more than two full- 
band counts (solid contours; see §4.3 for details). Dashed contours indicate 
the results obtained when only the 11 SDSS RQQs at z > 4.8 are considered 
in the joint spectral fitting. 

On average, it appears that any neutral X-ray absorption 
in z > 4.8 X-ray detected SDSS quasars has a column den- 
sity of 7Yh < 3.95 x 10 22 cirr 2 (at the 90% confidence level; 

2 when only the RQQs are used in 
For comparison, X-ray absorption of 
2 has been detected in a few RQQs and 
■ 1-3 (e.g., Cappi et al. 1997; Reeves & 



N H < 4.37 x 10 zz cm 
the spectral analysis). 



(2-5) xlO 22 cm 



some RLQs at z 
Turner 2000). Given the limited counting statistics available 
for the present sample, it is not possible to obtain useful spectral 
constraints for more complex absorption models (e.g., ionized 
absorption). The typical upper limits on the equivalent widths 
of narrow iron Ka lines (either neutral or ionized) are in the 
range 420-950 eV. 

The X-ray spectral results at z > 4 presented here and in 
V03 allow us to place constraints on any spectral changes in 
the X-ray continua of RQQs as a function of redshift. In Fig- 
ure 8 we plot the photon index as a function of redshift for some 
optically selected samples of RQQs. The photon indices plot- 
ted in the figure are for the nuclear, high-energy ( £ 2 keV) 
component, in order to avoid spectral complexities such as 
soft excesses. The ASCA RQQs have been taken from Vig- 
nali et al. (1999; V99 in the key), George et al. (2000; GOO), 
Reeves & Turner (2000; RT00), and Vignali et al. (2001b; 
VOlb). The BeppoSAX RQQs are from Mineo et al. (2000; 
MOO). At the highest redshifts we plot the z = 3.91 gravitation- 
ally lensed BALQSO APM 08279+5255 (Chartas et al. 2002; 
C02), the results obtained from joint X-ray spectral fitting of 
nine PSS quasars in the redshift range z ~ 4.1—4.5 (V03), and 
the 1 1 SDSS RQQs at z > 4.8 presented in this paper. The 
distribution of X-ray photon indices is characterized by signif- 
icant scatter, so precise parametric modeling of any F versus 
redshift relation is difficult. However, according to a Spear- 
man's p correlation test, a correlation between F and redshift 
is not significant (2.3a). Applying a linear regression algo- 
rithm that accounts for intrinsic scatter and measurement errors 
(Akritas & Bershady 1996), we obtain | |j | < 0.05 ±0.03. 
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FIG. 8. — Plot of photon index versus redshift for optically selected RQQs. 
The ASCA RQQs have been taken from Vignali et al. (1999; V99 in the key), 
George et al. (2000; GOO), Reeves & Turner (2000; RT00), and Vignali et al. 
(2001b; VOlb). The BeppoSAX RQQs are from Mineo et al. (2000; MOO). 
At the highest redshifts we plot the z = 3.91 gravitationally lensed BALQSO 
APM 08279+5255 (Chartas et al. 2002; C02), the results obtained from joint 
X-ray spectral fitting of nine PSS quasars in the redshift range z ~ 4.1-4.5 
(V03), and the 11 SDSS RQQs at z > 4.8 presented in this paper (the latter 
two data points are plotted at the average redshifts of the relevant papers). Er- 
ror bars are at the 68% confidence level for all data points. The object with the 
steepest X-ray photon index in this plot, PG 1543+489 (George et al. 2000), 
is known to be a narrow-line type 1 quasar (Boroson & Green 1992). 



5. SUMMARY 

We have reported the Chandra and XMM-Newton observa- 
tions of a sample of 13 quasars at z ~ 4.7-5.4, the majority of 
them selected by the SDSS, and compared their properties with 
those of all the z > 4 quasars previously observed in the X-ray 
band. The principal results are the following: 

• Eleven quasars have been detected in the X-ray band, 
doubling the number of z > 4.8 X-ray detected quasars. 
The X-ray detections include SDSS 0913+5919, the 
highest redshift (z = 5.11) RLQ thus far detected in 
the X-ray band. Two quasars, SDSS 0338+0021 and 
SDSS 1737+5828, were not detected by the present 
observations. The former non-detection is probably 
due to the short exposure time of the XMM-Newton 
observation after removal of flaring-background inter- 
vals. For the latter quasar, strong X-ray absorption 
(N H £ 1-2 x 10 23 cm" 2 at the source redshift) is a likely 
explanation. An HET spectrum of this source, revealing 
absorption features, seems to support this hypothesis. 
For our sample, no evidence for source extension or an 
overdensity of companion objects has been found. 

• The presence of a significant correlation between the 
soft X-ray flux and AZ?i45Q(i +z ) magnitude has been 
extended to the highest redshifts. This indicates that 
the same engine (namely accretion onto a supermassive 
black hole) is powering the ultraviolet and X-ray emis- 
sion. 
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• The anti-correlation between a ox and rest-frame 2500 A 
luminosity density found for SDSS RQQs over the red- 
shift range z ~ 0.2-6.3 (VBS03) has been confirmed 
(3.5-4.0a significance level) after the inclusion of the 
RQQs presented in this paper. Similarly to VBS03, 
no significant correlation between a ox and redshift has 
been found. These results suggest that the mechanism 
driving quasar broad-band emission is similar in the lo- 
cal and early Universe, with no evidence for unusual 
phenomena in high-redshift quasar nuclei such as "trap- 
ping radius" effects or accretion-disk instabilities. 

• The joint w 2.5-36 keV rest-frame X-ray spectrum of 
the 13 SDSS quasars at z > 4.8 observed thus far by 
Chandra is well parameterized by a simple power-law 
model with Y = 1.84+Jj ^ and no evidence for intrinsic 
absorption (N H ~ 4.0 x 10 22 cm~ 2 at 90% confidence). 

• There is no evidence for significant (> 30%) optical 
variability over the time interval of a few years (in the 
observed frame) between the SDSS and HET observa- 
tions. 

Among the quasars presented in this paper, the only 
object suitable for moderate-quality (rs 1000 photons) 
X-ray spectroscopy in a 50 ks XMM-Newton observation 
is SDSS 0231-0728. X-ray spectroscopy of the other 
quasars will be possible only with future missions such as 
Constellation-X, XEUS, and Generation-X (see V03 for some 
relevant simulations). 
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APPENDIX 

Here we report the updated X-ray flux values (and derived quantities) for z > 4 quasars observed by Chandra in the past. The 
X-ray values have been corrected for the ACIS quantum energy decay at low energies. 



Table A 1 

Optical, X-ray, and Radio Properties of z > 4 Quasars Previously Observed by Chandra 



Object 


z 


AB 1450(l+z) 


■^2500 A 




M B 


A.5-2 kcV 


fl keV 


log(L 2 kcv) b 


lOg(L 2 _10kev) 




R 


PSS 0059+0003 


4.16 


19.3 


1.06 


31.54 


-27.62 


10.2 


7.88 


27.41 


45.30 


-1.59 


<3.3 


BR 0103+0032 


4.43 


18.7 


1.85 


31.82 


-28.33 


20.1 


16.28 


27.77 


45.66 


-1.56 


2.9 


PSS 0121+0347 


4.13 


18.5 


2.22 


31.86 


-28.41 


41.3 


31.65 


28.01 


45.90 


-1.48 


300.3 


PSS 0133+0400 


4.15 


18.2 


2.93 


31.98 


-28.72 


20.8 


16.01 


27.72 


45.61 


-1.64 


<4.4 


PSS 0134+3307 


4.53 


18.5 


2.22 


31.92 


-28.57 


11.5 


9.45 


27.55 


45.44 


-1.68 


<6.1 


PSS 0209+0517 


4.14 


17.8 


4.24 


32.14 


-29.12 


15.2 


11.66 


27.58 


45.47 


-1.75 


<3.0 


SDSS 0210-0018 


4.77 


19.3 


1.06 


31.63 


-27.85 


37.3 


32.12 


28.11 


46.00 


-1.35 


86.1 


SDSS 0211-0009 


4.90 


20.0 


0.56 


31.37 


-27.20 


3.8 


3.36 


27.15 


45.04 


-1.62 


<7.5 


BR 0241-0146 


4.06 


18.4 


2.44 


31.88 


-28.48 


4.3 


3.26 


27.01 


44.84 


-1.87 


< 1.5 


PSS 0248+1802 


4.43 


18.1 


3.21 


32.06 


-28.93 


21.2 


17.16 


27.79 


45.68 


-1.64 


<4.2 


BR 0308-1734 


4.00 


18.1 


3.21 


31.99 


-28.76 


4.1 


3.06 


26.97 


44.86 


-1.93 


< 3.9 


BR 0401-1711 


4.23 


18.9 


1.54 


31.71 


-28.05 


12.9 


10.07 


27.53 


45.42 


-1.61 


< 8.5 


SDSS 0836+0054 


5.82 


18.8 


1.67 


31.96 


-28.67 


12.3 


12.50 


27.83 


45.73 


-1.58 


9.8 


PSS 0926+3055 


4.19 


16.7 


11.67 


32.59 


-30.24 


33.5 


25.94 


27.93 


45.82 


-1.79 


<0.3 


PSS 0955+5940 


4.34 


18.4 


2.44 


31.93 


-28.60 


6.0 


4.79 


27.22 


45.12 


-1.81 


<5.4 


PSS 0957+3308 


4.20 


18.2 


2.93 


31.99 


-28.74 


11.3 


8.78 


27.46 


45.36 


-1.74 


< 1.1 


SDSS 1030+0524 


6.28 


19.7 


0.76 


31.67 


-27.94 


2.5 


2.73 


27.22 


45.12 


-1.71 


<6.8 


BR 1033-0327 


4.51 


18.8 


1.69 


31.80 


-28.26 


8.4 


6.87 


27.41 


45.30 


-1.69 


<2.5 


PSS 1057+4555 


4.12 


17.6 


5.09 


32.22 


-29.31 


25.3 


19.34 


27.80 


45.68 


-1.70 


2.3 


SDSS 1129-0142 


4.85 


19.2 


1.17 


31.68 


-27.98 


<2.4 


< 1.75 


< 26.86 


< 44.84 


< -1.85 


<3.4 


SDSS 1204-0021 


5.03 


19.1 


1.28 


31.75 


-28.14 


4.6 


4.13 


27.26 


45.15 


-1.72 


< 3.4 


SDSS 1208+0010 


5.27 


20.5 


0.35 


31.22 


-26.82 


1.3 


1.24 


26.77 


44.66 


-1.71 


< 13.1 


PC 1247+3406 


4.90 


19.2 


1.17 


31.69 


-28.00 


11.0 


9.66 


27.61 


45.50 


-1.57 


< 3.3 


SDSS 1306+0356 


5.99 


19.6 


0.85 


31.68 


-27.97 


5.7 


5.94 


27.53 


45.42 


-1.59 


<5.3 


PSS 1317+3531 


4.36 


18.9 


1.54 


31.73 


-28.10 


3.6 


2.91 


27.01 


44.90 


-1.81 


<2.3 


PSS 1326+0743 


4.09 


17.2 


7.36 


32.37 


-29.70 


29.4 


22.32 


27.85 


45.74 


-1.73 


<0.5 


PSS 1347+4956 


4.51 


17.4 


6.12 


32.36 


-29.66 


17.5 


14.35 


27.73 


45.62 


-1.78 


0.1 


PSS 1435+3057 


4.35 


19.1 


1.28 


31.65 


-27.90 


<3.6 


<2.87 


< 27.00 


< 44.90 


< -1.79 


<2.9 


PSS 1443+2724 


4.42 


19.0 


1.40 


31.70 


-28.03 


11.3 


9.10 


27.51 


45.41 


-1.61 


<2.8 


PSS 1443+5856 


4.26 


17.7 


4.64 


32.20 


-29.26 


4.7 


3.66 


27.09 


44.99 


-1.96 


<2.8 


SDSS 1532-0039 


4.62 


19.4 


0.97 


31.57 


-27.70 


<2.5 


<2.10 


< 26.91 


< 44.80 


< -1.79 


<4.3 


SDSS 1605-0112 


4.92 


19.4 


0.97 


31.61 


-27.80 


<3.0 


<2.62 


< 27.05 


< 44.94 


< -1.75 


<4.5 


BR 2212-1626 


4.00 


18.6 


2.03 


31.79 


-28.25 


13.4 


9.95 


27.48 


45.38 


-1.65 


<6.2 



Note. — Units are the same as in Table 3. 
"Rest-frame 2500 A luminosity density (erg s _1 Hz -1 ). 
b Rest-frame 2 keV luminosity density (erg s _1 Hz -1 ). 



